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PREFACE 


The  research  reported  herein  was  concl\  cted  by  the  Geophysics 
Group  of  the  Infrared  and  Optics  Division  of  the  Environmental  Research 
Institute  of  Michigan. 

This  work  was  sponsored  by  the  Naval  Regional  Procurement  Office 
Long  Beach,  California  under  Contract  N00123-74-C-0761  and  monitored  by 
the  Naval  Electronic  Laboratory  Center  (NELC) . The  period  of  research 
extended  from  26  November  1973  to  15  December  1975.  Project  monitor 
on  this  contract  was  John  G.  Hoffman  (NELC). 

The  principal  investigator  was  Fred  J.  Tanis  with  important 
contributions  to  the  technical  program  made  by  J. Adams.  This  research 
was  guided  by  Mr.  R.R.  Legault,  Director  of  the  Infrared  and  Optics 
Division,  The  Institute  report  number  i§  IO68OO70-F. 
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SUMMARY 


This  contract  was  directed  toward  development  of  a seismic  mortar 
location  system.  The  concept  envisions  a small,  low-cost,  man  packed 
computer,  and  hand-implaced  sensors  with  data  telemetry.  The  system 
operator  acts  as  a forward  observer  calling  fire  on  enemy  gun  positions. 

The  sensor  packages  would  have  threshold  detection,  pass  band  filtering, 

A/D  conversion,  discrimination  logic,  and  data  transmission  capabilities. 
Target  positioning  calculations  and  signal  analysis  would  be  performed 
in  a small  computer.  Operator  input  includes  sensor  and  calibration  shot 
coordinates.  Output  consists  of  status  indicators  and  enemy  gun  location 
estimates  which  can  be  transmitted  to  fire  control.  This  system  would 
also  be  applied  to  locate  large  artillery  and  other  targets  with  impulsive 
seismic  signatures. 

The  principal  advantages  of  a seismic  system  over  the  sound  ranging 
approach  are  (1)  the  use  of  velocities  which  are  not  time  dependent 
allowing  adaptive  location,  (2)  terrain  features  will  not  mask  the  seismic 
signals  as  with  acoustics,  and  (3)  seismic  signals  allow  direct  discrimina- 
tion of  multiple  events  and  signal  identification  through  match  comparisons. 
The  principal  disadvantages  of  a seismic  locater  as  with  all  passive 
systems  is  a noise  limited  range.  This  study  found  that  further  effort 
must  be  directed  towards  signal  enhancement. 

Seismic  surface  waves  constitute  the  largest  signal  at  the  ranges 
under  consideration  and  were,  therefore,  used  as  a basis  for  developing 
a location  technique. 

Two  location  algorithms  were  derived  to  utilize  seismic  data.  The 
first  was  similar  to  hyperbolic  location  but  used  a Kalman  filter  to 
correct  for  individual  path  velocities.  The  second  technique  compared 
signals  from  enemy  recoil  with  those  of  directed  return  fire  in  a "seismic 
on  seismic"  location  estimate.  Computer  simulation  of  derived  location 
algorithms  was  found  useful  to  obtain  capabilities  in  terms  of  error 
characteristics . 
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A field  measurement  program  was  initiated  to  provide  a data  base 
for  further  evaluation.  Seismic  data  recordings  were  made  of  mortar 
recoils  and  explosives  at  Twenty  Nine  Palms  and  Camp  Pendleton  Marine 
bases. 

The  Camp  Pendleton  data  were  of  a poor  quality  because  of  local 
noise  problems  and  as  such  did  not  produce  useful  results  when  applied 
to  the  location  problem.  The  TWenty  Nine  Palms  experiment  on  the  other 
hand  produced  good  data.  While  observations  were  made  at  ranges  less 
than  2.0  kilometers,  it  is  felt  that  seismic  recoil  signals  without 
enhancement  would  be  useful  under  quiet,  conditions  to  at  least  3.0 
kilometers.  Signal  enhancement  processing  could  double  the  effective 
range. 

Location  accuracies  produced  by  windless  acoustic  data  (an  ideal 
condition)  were  generally  better  than  those  from  corresponding  seismic 
data.  Because  of  the  higher  seismic  velocities  seismic  results  were 
more  sensitive  to  timing  errors.  Significant  advances  must  be  made  in 
the  methods  of  observing  seismic  wave  time  delays  before  the  advantages 
of  a seismic  location  system  can  be  realized.  Adaptive  processing  of 
seismic  signals  could  also  reduce  noise  without  any  intervention  by  the 
operator. 
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1 

INTRODUCTION 

The  problem  of  enemy  mortar  location  hap  not  been  solved  for  a 
small  forward  base  of  company  size.  The  cost  and  logistics  of  using 
radar  for  this  purpose  are,  to  date,  prohibitively  expensive.  This 
exploratory  work  as  part  of  the  U.S.  Marine  Corps  HOWLS  (Hostlle-Weapons- 
Location-System)  program  has  been  directed  toward  development  of  a 
seismic  mortar  location  system.  This  concept  envisions  development  of 
a small,  low  cost,  man-packed  computer,  and  hand-implaced  sensors  with 
data  telemetry  capability.  The  system  operator  will  function  as  a forward 
observer  calling  fire  for  calibration  and  on  enemy  gun  positions. 

Unlike  the  radar  in-flight  projectile  approach,  this  program  concentrates 
on  the  propagation  of  seismic  energy  resulting  from  the  impulse  of  the 
firing  weapon  and  could  provide  360°  surveillance  cover.  The  sensor 
packages  which  would  probably  be  expendable  are  envisioned  to  have 
threshold  detection,  passband  filtering,  A/D  conversion,  a small 
amount  of  digital  storage,  limited  discrimination  logic,  and  data 
transmission  capability.  All  target  positioning  calculation  and 
signal  analysis  will  be  performed  in  a small  computer.  Operator  input 
will  include  sensor  and  calibration  shot  position  information.  Output 
will  consist  of  status  indicator  and  enemy  position  information  which 
can  be  transmitted  to  fire  control. 

One  of  the  principal  advantages  of  the  seismic  system  over  the 
conventional  sound  ranging  system  is  the  use  of  seismic  velocities 
which  are  not  time  dependent.  This  feature  allows  the  location  to  take 
on  an  adaptive  character.  With  repetition  of  recoils  and  impacts  more 
knowledge  of  media  path  velocities  is  obtained  with  which  the  position 
of  enemy  fire  can  be  progressively  estimated.  Because  a seismometer 
has  some  capability  as  a microphone,  acoustic  muzzle  blast  signals  can 
also  be  utilized,  if  needed,  for  location. 
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Terrain  features  wi3.1  not  mask  the  seismic  signal  from  a muzzle 
blast  as  they  will  for  acoustics.  As  return  fire  is  ranged  upon  the 
enemy's  gun  the  signatures  from  friendly  impact  and  enemy  recoil  will 
look  similar  and  can  be  incorporated  into  a "seismic  on  seismic" 
l.ocatlon  estimate.  Thus  we  will  know  convincingly  when  the  enemy 
fire  has  been  neutralized.  The  practical  use  of  these  techniques  in 
a mortar  location  algorithm  has  been  investigated. 

Discrimination  of  seismic  recoil  signals  from  shell  impact 
signals  is  necessary  when  multiple  events  occur  with  rapid  succession. 
The  seismic  mortar  recoil  signals  are  repeatable  which  allows  signal 
identification  through  successive  match  comparisons.  In  addition, 
directionality  of  the  acoustic  wave  could  possibly  discri.ninate 
signals  from  Improbable  azimuths. 

As  with  all  passive  systems  a seismic  mortar  locator  will  be 
noise  limited.  If  noise  conditions  are  very  local  then  perhaps  the 
problem  can  be  relieved  by  ignoring  one  or  more  sensors.  A generally 
high  battlefield  noise  condition  will  reduce  the  effective  range. 

The  approach  taken  to  investigate  the  feasibility  of  a seismic 
mortar  location  sj'ctem  included  the  following  program  steps. 

1.  Utilizing  existing  data  examine  possible  mortar  generated 
signals  and  characteristics  which  could  be  exploited  in  a 
mortar  location  scheme  with  attention  to  problems  of  detec- 
tion, discrimination,  and  location. 

2.  Develop  a mortar  location  algorithm  which  could  be  used  in  a 
light  weight  man-packed  system.  Simulate  location  algorithms 
and  evaluate  their  expected  performance  on  assumed  input 
errors . 

3.  Conduct  a field  measurements  program  to  collect  data  which 
can  be  used  to  evaluate  the  algorithms. 

4.  Process  the  seismic  field  data,  analyze  its  characteristics 

and  utilize  it  as  much  as  possible  to  obtain  locations  of 
mortar  fire. 
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SEISMIC  MORTAR  SIGNATURE  CHARACTERISTICS 

A mortar  firing  generates  four  signals  which  will  be  recorded 
by  the  seismometer.  These  include  seismic  waves  generated  by  recoil, 
the  air  wave  from  the  muzzle  blast,  the  seismic  waves  from  the 
explosive  shell  impact,  and  the  air  wave  from  the  shell  burst. 

A shell  burst  event  (an  explosion)  generates  seismic  energy  which 
penetrates  the  surface  layers  and  propagates  to  the  sensor  as  a 
compressional  wave  (P  wave)  or  simply  as  sound  in  earth.  Also 
generated,  and  perhaps,  with  the  greatest  amount  of  energy,  are  surface 
waves  (Rayleigh  wave) . From  such  an  event  very  little  seismic  shear 
wave  energy  is  observed.  In  contrast,  mortar  recoil  will  generate 
shear  energies  confined  to  Love  modes  (surface  confined  shear  waves) 
and  shear  body  modes  which  can  be  polarized  in  both  the  horizontal  (SH) 
and  vertical  (SV) . Observed  polarization  will  be  partly  determined 
by  azimuth  to  the  mortar  location. 

Seismic  signal  from  mortar  recoil  has  been  shown  to  be  highly 
repeatable.  Signals  received  from  explosions,  on  the  other  hand, 
while  showing  a great  deal  of  correlation,  do  not  show  the  same  degree 
of  repeatability  unless  they  are  small  explosions  and  detonated  at  the 
same  position.  It  is  not  the  kind  of  event  but  rather  its  degree  of 
simularity  which  determines  repeatability.  When  explosions  or  impacts 
are  scattered  about  due  to  range  and  deflection  errors  they  exhibit 
less  correlation.  Such  tendencies  seem  reasonable  since  shells  explode 
at  different  distances  and  azimuths  from  the  sensor. 

Such  a discriminant  could  be  implemented  by  storing  a known 
recoil  signal  from  enemy  mortar  along  with  its  auto-correlation 
function.  A second  signal  would  then  be  compared  to  the  first  by 
ratioing  the  cross-correlation  peaks  with  the  stored  auto-correlation 
function.  A constant  ratio  would  indicate  a recoil  from  the  same 
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source.  Conversely,  a variable  ratio  would  indicate  recoil  from 
another  source  or  shell  burst. 

There  is  very  little  vibual  difference  in  the  nature  of  the  air 
wave  train  as  recorded  from  a shell  burst  or  a muzzle  blast.  In 
addition,  the  acoustic  spectra  of  these  air  waves  are  much  the  same. 

Thus  using  only  the  acoustic  air  waves,  it  may  be  difficult  if  not 
impossible  to  directly  discriminate  a muzzle  blast  from  a shell  burst. 

In  this  case,  one  is  left  to  sorting  out  all  the  arrivals  in  a large 
computer  in  order  to  eliminate  all  false  locations.  On  the  other 
hand,  the  seismic  waves  from  an  explosion  and  a recoil  may  differ 
substantially  in  terms  of  signal  level  on  each  of  the  three-component 
seismometers . 

2.1  ANALYSIS  OF  THE  AIR  WAVE 

The  air  wave  portion  of  the  muzzle  blast  seismogram  will  be 
the  largest  signal  on  the  trace  and  could  be  very  useful  in  determining 
azimuth  to  the  mortar  location.  Two  modes  of  seismic  signal  are 
possible.  First,  the  sound  wave  may  couple  directly  into  the  seismo- 
meter housing  or  immediate  ground  area  as  structural  borne  sound. 

Second,  the  sound  may  couple  into  the  ground  continuously  as  it 
propagates  to  the  seismometer  reinforcing  a Rayleigh  wave  v/hich 
propagates  at  the  velocity  of  sound.  This  wave  is  often  called  the 
air-coupled  Rayleigh  wave.  The  recorded  signal  can  be  dominated  by 
either  one  of  these  air-coupling  phenomenon.  The  prominence  of  an  air- 
coupled  Rayleigh  wave  will  depend  on  whether  such  a wave  can  be 
supported  by  the  ground  media.  If  the  layered  structure  in  the  earth 
will  support  such  a wave  in  the  frequency  band  of  the  acoustic  air 
wave  then  the  air-coupled  Rayleigh  wave  will  dominate. 

The  two  horizontal  seismometers  of  a three-component  package  can 
be  used  to  obtain  direction  such  as  one  might  with  a pair  of  directional 
microphones,  i.e.,  by  the  ratio  of  signal  levels  or  by  timing  the  waves 
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to  separated  stations.  If  the  shape  of  the  acoustic  wavefront 
becomes  distorted  by  the  presence  of  foliage,  and  atmospheric 
conditions,  good  azimuth  resolution  may  be  impossible.  The  seismic, 
direct-coupled  acoustic  wave  would  indicate  a comparable  degree  of 
distortion.  Thus  in  the  case  of  wave  distortion,  using  the  seismometer 
as  a micrphone  will  not  improve  azimuth  resolution.  The  air-coupled 
Rayleigh  wave  will  be  less  likely  distorted  because  its  energy  is 
coupled  to  the  ground  over  the  entire  wave  path.  Azimuthal  information 
from  the  acoustic  wave  could  be  used  as  (1)  an  initial  rough  location 
estimate,  and  ^2)  an  event  discriminant  for  unlikely  azimuths. 


2.2  DIRECTION  FINDING  WITH  ACOUSTIC  WAVES 

The  Investigation  has  been  made  of  using  two  horizontal  seismometers 
to  obtain  direction,  such  as  might  be  done  with  directional  micro- 
phones. A general  direction  can  be  found  from  the  acoustic  waves 
accompanying  a muzzle  blast.  Further  analysis  could  be  performed  to 
refine  angle  estimates  and  identify  major  causes  of  error.  It  has 
not  been  ascertained,  for  example,  whether  ripples  in  the  sound  front 
in  the  atmosphere  appear  in  the  Rayleigh  wave  front  in  the  earth. 

If  these  ripples  do  not  appear  in  the  ground  then  a major  source  of 
acoustic  direction  finding  error  will  be  eliminated  [1].  Analyses 
which  have  been  completed  for  this  study  are  reported  in  Section  4. 


2.3  USE  OF  SEISMIC  SURFACE  WAVES 

The  seismic  surface  waves  generated  by  recoils  and  impacts  may 
be  the  largest  signal  on  the  seismogram  outside  of  the  air  wave.  Body 
waves  will  generally  attenuate  rapidly  because  of  their  high  frequency 
character  and  spherical  spreading. 

As  with  the  acoustic  waves  these  surface  waves  may  be  used  to 
locate  enemy  fire  by  observing  the  arrival  time  difference  (AT's)  at 
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three  or  more  stations.  However,  accurate  hyperbolic  fixes  will 
require  knowledge  of  the  path  velocities.  The  Rayleigh  and  Love 
wave  velocities  can  be  expected  to  change,  not  only  with  area  of 
deployment,  but  also  to  a lesser  extent  within  the  actual  deployment 
area  itself.  The  approximate  range  and  azimuth  of  all  return  rounds 
will  be  known  and  the  seismic  waves  generated  by  these  impacts  and 
recoils  can  be  used  to  more  accurately  re-estimate  the  surface 
velocities.  Thus  a seismic  mortar  location  system,  once  deployed, 
will  improve  in  effectiveness  as  the  position  is  defended. 

2.4  DELAY  TIME  ANALYSIS 

Efforts  were  made  to  evaluate  some  existing  seismic  mortar  data 
for  the  use  of  correlation  for  AT  determinations.  Analog  mortar  and 
explosion  data  used  were  recorded  at  Ft.  Sill,  Oklahoma  under  Army 
contract  DAAK02-73-C-0230.  Selected  records  were  filtered  and 
digitized.  A simple  correlation  program  was  written  to  perform  the 
necessary  computations. 

Figure  1 shows  a comparison  of  three-component  records  of  an 
explosion  recorded  at  two  0.5  km  sites  which  are  separated  by  382.7  m. 
The  peak  of  each  correlation  function  was  determined  so  that  the  lag 
time  from  zero  time  alignment  could  be  determined.  This  quantity  has 
been  labeled  peak  error  and  assumes  path  velocities  are  Identical. 

In  the  figure  Rayleigh  waves  which  appear  on  the  vertical  and  radial 
channels  appear  to  have  similar  path  effects.  The  Love  waves  which 
appear  on  the  transverse  channel  correlate  poorly  and  produce  a 
large  alignment  error.  Additional  records  examined  at  the  same  site 
indicated  that  Rayleigh  waves  which  had  travelled  along  substantially 
different  paths  remained  well  correlated.  The  results  from  another 
site  at  Ft.  Sill  were  not  as  encouraging.  Figures  2 and  3 show  the 
same  type  of  comparison  for  a mortar  signature.  The  results  were 
generally  unacceptable  which  may  be  due  to  the  presence  of  special 
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Vertical  Component 

SEISMICS  ACOUSTICS 


Correlation  Function  Peak  Error  87  ms 


Transverse  Component 

SEISMICS  ACOUSTICS 


Site  4 
Site  3 


Correlation  Function  Peak  Error  9ms 


Radial  Component 


FIGURE  1.  COMPARISON  OF  THREE  COMPONENT  SEISMIC  SIGNATURES 
FROM  A TNT  SHOT  AT  TWO  SITES  SEPARATED  BY 
382,7  m. , Fort  Sill,  Oklahoma,  1971 
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Correlation  Function  wiKlIlln  Peak  Error  18ms 
Vertical  Component 


Site  0 
Site  1 


Correlation  Function  Peak  Error  90ms 

Transverse  Component 
SEISMICS 


ACOUSTICS 


Correlation  Function  Peak  Error  97ms 

Radial  Component 

FIGURE  2.  COMPARISON  OF  THREE  COMPONENT  SEISMIC  MORTAR  SIGNATURES 
SEPARATED  BY  382.7  m. , Fort  Sill,  Oklahoma,  1971 
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SEISMICS 


Correlation  Function  Peak  Error  27ms 

Vertical  Component 


Site  7 
Site  0 

Correlation  Function  Peak  Error  66ms 

Transverse  Component 

SEISMICS  ACOUSTICS 


Correlation  Function  • • Peak  Error  57ms 

Radial  Component 

FIGURE  3.  COMPARISON  OF  THREE  COMPONENT  SEISMIC  MORTAR  SIGNATURES 
SEPARATED  BY  382.7  m.  , Fort  Sill,  Oklahoma,  1971 
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geologic  features.  Seismic  record  comparison  of  a mortar  recoil  with 
an  explosion  detonated  at  the  firing  position  are  shown  in  Figure  4. 
Again,  excellent  correlation  alignment  was  achieved  on  the  vertical 
and  radial  channels  (due  to  Rayleigh  waves)  but  remained  poorly 
correlated  on  transverse  channels.  In  addition,  it  is  evident  that 
the  recoil  generates  much  larger  Love  waves  than  from  the  explosion. 
This  large  difference  may  account  for  the  large  peak  lag  error. 


2.5  SENSOR  REQUIREMENT 

In  considering  deployment  of  the  seismic  mortar  location  system 

for  point  defense  one  may  ask  how  many  sensors  are  required  for  360° 

coverage.  If  each  sensor  has  a reliable  detection  range  of  r and 

max 

it  is  desired  to  locate  enemy  fire  from  a radius  of  R around  the  point 
defense  than  the  sensor  packages  would  optimumly  be  placed  uniformly 
on  a circle  of  some  radius  a.  Sensors  on  the  opposite  side  from 
the  enemy  mortar  will  be  out  of  detection  range  for  radii  greater  than 
^max”^'  Suppose,  further,  that  we  require  a minimum  of  k sensors  to 
detect  the  enemy  mortar  recoil.  Since  the  k sensors  will  all  lie  on 
a circular  arc,  the  requirement  reduces  to  having  the  two  end  positions 
at  maximum  range.  Of  the  arc  subtends  an  angle  iji  then  N > ZTtk/iJi 
sensors  are  required.  The  relationship  between  N and  is  easily 

derived  as 


N = Tr(k-l)/arc  cos  [ (R^-  + - r^  )/2  Ra] 

max 


Under  this  formulation  N is  very  sensitive  to  r .A  more  useful 

max 

relation  is 


r = / (R  - a cos 


i + a2.i„2  lltli 


max 


N 


N 
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Correlation  Function  a0|^VA  Peak  Error  Sms 

Vertical  Component 


Correlation  Function  * Peak  Error  66ms 

Transverse  Component 


FIGURE  4.  COMPARISON  OF  THREE  COMPONENT  SEISMIC  SIGNATURES  FROM 

TNT  SHOT  WITH  A MORTAR  RECOIL.  Fort  Sill,  Oklahoma,  1971 
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If  unlimited  numbers  uf  sensors  were  available  then  by  the  above 
equation 


Thus  we  can  examine  N such  that  r begins  to  deviate  rapidly  from 

max 

R-a.  Using  a value  of  0.5  and  1.0  km,  an  R value  of  3.0  km,  and  k = 3 
the  point  of  ten  percent  departure  from  R -a  ■ 2.0  lies  between 
8 and  10  sensors.  For  a k value  of  4 these  numbers  Increase  to  12 
and  15  sensors  respectively. 

Examination  of  this  circular  criterion  for  sensor  placement 
leads  under  practical  considerations  to  a conclusion  that  a minimum 
of  10  sensors  are  needed  for  360“  coverage. 

2.6  SIGNATURE  DISCRIMINATION  PROBLEM 

As  indicated  previously  discrimination  of  seismic  recoil  signals 
from  impact  signals  may  be  crucial  to  effective  operation  of  a mortar 
location  system.  A discrimination  capability  is  especially  necessary 
during  periods  of  rapid  exchange  of  fire.  Lack  of  such  capability 
necessitates  a lengthy  and,  perhaps,  complex  sorting  problem. 

Examination  of  seismic  data  collected  at  Ft.  Sill  during 
November  1973  and  under  this  project  suggest  three  possible  approaches 
for  seismic  discrimination  of  mortar  recoil  and  Impact  signatures. 

Results  from  this  study  have  suggested  three  potential  discrimi- 
nants . 

1.  Seismic  signal  for  mortar  recoil  has  been  shown  to  be 
highly  repeatable.  Signals  received  from  explosions,  on  the  other 
hand,  while  showing  a great  deal  of  correlation,  do  nor  show  the 
same  degree  of  repeatability  unless  they  are  small  explosions  and 
detonated  at  one  position.  When  explosions  or  impacts  are  scattered 
about  due  to  range  and  deflection  errors  they  exhibit  less  correlation. 


12 


rn 


Terim  

FORMiSl^VILLOW  run  LAaORATORIKS^THC  UNiVCntn77^!MCHloAN 


Such  tendencies  seem  reasonable  since  shells  explode  at  different 
distances  and  azimuths  from  the  sensor.  Figures  5 and  6 show  a 105  mm 
shell  impact  signature  as  recorded  at  0.5  km  distance.  While  some  of 
the  surface  waves  show  correlation  the  general  repeatability  of  the 
record  is  much  less  than  in  the  case  of  recoil.  Thus  repeatability 
of  the  signal  shows  promise  as  a signal  discriminant. 

2.  Explosive  sources  tend  to  produce  relatively  stronger  body 
wave  energy  than  recoil  sources.  A body  wave/surface  ratio  is  then 

a potential  discriminant  between  these  two  sources.  The  discriminant 
ratio  may  be  enhanced  by  comparing  signals  from  selected  narrow  band 
filters.  While  this  discriminant  may  work  well  for  close  distances 
where  body  wave  signals  from  mortar  can  be  observed,  there  are  obvious 
difficulties  at  the  far  distances  where  they  are  highly  attenuated 
and  difficult  to  measure. 

3.  Since  surface  waves  attenuate  less  rapidly  than  body  waves 
they  can  be  observed  at  greater  distances.  Thus  a surface  wave 
discriminant  would  have  more  potential  use  for  the  mortar  location 
problem.  A procedure  was  developed  to  look  at  the  sum  of  squares 

of  the  longitudinal  and  vertical  components.  Strong  phase  correlation 
was  observed  for  surface  waves  recorded  beyond  1 km  from  explosions 
at  about  25  Hz.  For  surface  waves  from  mortar  recoil  phase  correla- 
tion was  observed  to  vary  little  fiom  2 Hz  to  50  Hz.  The  presence  or 
absence  of  such  phase  correlation  could  potentially  be  used  to 
discriminate  recoils  from  impacts. 

It  was  first  thought  that  recoil  seismics  produced  less  energy 
on  the  vertical  component  than  explosions.  If  this  were  the  case,  an 
energy  ratio  between  components  could  be  used  as  a discriminant. 
However,  when  tested  this  discriminant  was  found  to  be  unreliable. 

Of  the  discriminants  considered  repeatability  appears  to  be  the  most 
practical  approach.  Such  a discriminant  could  be  implemented  by 
storing  a known  recoil  signal  from  enemy  mortar  along  withs  its 
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Figure  5.  Three  Component  Seismic  Signatures, 
105  mm  Shell  Impact* Recorded  at 
Ft.  Sill,  11/11/73,  Site  4 @ 736. 
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auto-correlation  function.  A second  signal  would  then  be  compared 
to  the  first  by  ratioing  the  cross-correlation  peaks  with  the  stored 
auto-correlation  function.  A constant  ratio  would  indicate  a recoil 
from  the  same  source.  Conversely,  a variable  ratio  would  indicate 
recoil  from  another  source  or  shell  burst. 

2.7  DISPERSION  OF  SURFACE  WAVES 

The  Rayleig.;  and  Love  surface  modes  generated  by  mortar  recoil 
and  shell  impacts  will  travel  at  velocities  which  vary  with  frequency. 

If  the  media  is  dispersive  such  velocity  differences  will  be  great 
and  the  surface  wave  train  will  appear  spread  out.  Such  phenomena 
could  greatly  effect  the  accuracy  of  time  delay  correlations  if  the 
path  lengths  or  media  characteristics  are  different.  Early  in  this 
program  model  studies  using  multi-layered  dispersion  equations  were 
initiated  to  evaluate  the  problem  of  dispersion. 

Mortar  recoil  tends  to  generate  Love  surface  mode  while  shell 
impacts  generate  predominantly  Rayleigh  mode.  If  arrivals  from  impacts 
are  to  be  used  to  correct  times  for  mortar  recoil  or  if  they  are  to 
be  incorporated  into  a seismic  on  seismic  correction  of  return  fire 
then  such  velocities  must  be  similar  and  in  constant  ratio  over  the 
principal  frequency  components.  Data  recorded  for  the  field  tests 
conducted  at  Ft.  Sill,  Oklahoma,  October  1973  wave  has  been  analyzed  for 
consistency  of  explosions  and  recoils.  In  addition,  using  a 
refraction  velocity  profile  obtained  from  the  Ft.  Sill,  1971  data 
theoretical  models  of  surface  wave  velocity  were  calculated  for 
various  modes. 

A series  of  equations  [2],  describing  the  dispersion  of  surface 
waves  on  multi-layered  media,  were  programmed  for  the  digital 
computer.  They  relate  surface  wave  phase  velocity  to  frequency. 

The  program  input  consists  of  a layered  model  of  shear  (3)  and 
compressional  (a)  wave  velocities  along  with  density  ratios. 
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Figures  7 and  8 show  the  resulting  phase  velocity  curves  for  both  Love 
and  Rayleigh  waves.  The  first  four  modes  have  been  determined.  The 
curves  have  roughly  the  shape  of  a hyperbola  except  that  the  Love 
modes  have  as  an  upper  limit  to  their  phase  velocities  which  is  the 
shear  velocity  in  the  second  layer  or  half-space.  The  compressional 
velocities  are  based  on  velocity  data  obtained  at  Ft.  Sill  and  shear 
velocities  were  obtained  from  the  Handbook  [3].  The  only  difference 
in  Models  I and  II  is  that  Model  II  has  a layer  of  twice  the  thickness 
of  Model  I.  The  effect  of  thickening  the  layer  is  to  lower  the 
frequency  at  which  there  is  a strong  upward  trend  of  the  phase 
velocity  curve.  The  upper  and  lower  limits  are  controlled  by  the 
layer  and  half  space  shear  velocities  respectively. 

In  general,  the  models  indicate  that  Love  modes  are  close 
approximations  to  the  Rayleigh  modes  for  high  frequencies,  i.e., 
the  flat  portions  of  the  various  curves.  The  worst  approximations 
occur  below  frequencies  where  the  mode  curves  turn  upward.  This 
frequency  corresponds  to  a wavelength  roughly  equal  to  one-fourth  the 
layer  thickness.  Thus  is  the  geology  has  significant  layering  with 
thickness  three  or  four  times  the  longest  wavelength  of  interest  we 
can  be  assured  that  Love  and  Rayleigh  velocities  are  nearly  the 
same.  For  lower  frequencies  Love  and  Rayleigh  velocities  will  have 
an  increasing  ratio. 

The  phase  velocity  curves  cannot  be  applied  directly  to  the 
available  data  because  shear  velocities  were  unknown.  In  addition 
higher  frequency  waves  tend  to  shift  predominant  energies  to  higher 
modes.  Mode  identification  is  not  possible  with  the  limited  Ft.  Sill 
data. 

The  use  of  impact  velocities  for  prediction  of  recoil  surface 
wave  velocities  appears  to  be  substantially  enhanced  if  we  obtain 
such  velocities  from  the  vertical  component.  Ft.  Sill  data  shows  a 
very  repeatable  vertical  component  from  small  explosions  to  recoils 
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Figure  7.  Theoretical  Rayleigh  and  Love  Surface  Wave  Dispersion 
Curves  for  Layer /Half-Space  Model  I. 


Theoretical  Rayleigh  and  Love  Surface  Wave  Dispersion 
Curves  for  Layer/Half-Space  Model  II. 
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over  the  same  path  (Figure  A).  As  expected  the  transverse  channels 
show  the  greatest  dissimilarity  between  explosions  and  recoils.  The 
explanation  may  be  that  Rayleigh  and  Love  modes  are  generated  In 
some  degree  by  both  sources  and  that  Rayleigh  modes  show  strongest 
motions  on  the  vertical  and  longitudinal  channels  while  Love  waves 
exhibit  little  motion,  if  any,  on  the  vertical  channel. 

The  model  studies  which  were  made  indicate  that  dispersion  effects 
are  not  significant  over  the  short  distance  changes  to  be  experienced 
in  impact/recoil  analysis.  Dispersion  effects  can,  however,  be 
significant  when  comparing  signals  which  have  travelled  over  great 
differences  in  distance.  It  should  be  possible  to  utilize  this 
characteristics  for  discrimination  or  cataloging  of  events  which 
occur  over  a broad  range  of  distances  from  the  recording  sensor 
package . 
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DEVELOPMENT  OF  A LOCATION  ALGORITHM 

The  location  algori-hm  consists  of  those  processing  elements  and 
mathematical  computing  steps  needed  to  obtain  a mortar  location  esti- 
mate from  sensor  time  of  arrival  delays.  In  this  program  two 
computational  algorithms  for  location  were  developed.  Both  techniques 
were  based  on  the  observation  that  mortar  recoils  and  impacts  both 
produce  simular  seismic  signals  as  measured  by  a vertical  axis 
seismometer.  However,  these  event  types  do  look  appreciably  different 
when  observed  in  the  horizontal  plane  due  to  the  fact  that  gun 
recoil  generates  substantial  Love,  (shear)  surface  waves  and  an  impact 
explosion  generates  only  very  limited  amounts.  The  vertical  axis 
signals  are.  dominated  by  Rayleigh  surface  wave  energy.  Because 
these  waves  have  the  largest  amplitudes  they  can  usually  be  observed 
at  greater  distances  from  equal  noise  levels  then  first  arrival 
compressional  waves.  The  exception  to  this  case  would  occur  if  the 
media  were  sufficiently  dispersive  so  as  to  spread  out  the  surface 
wave  train  causing  a rapid  decrease  in  amplitude. 

The  first  technique  described  below  does  not  rely  specifically 
on  either  the  surface  or  the  compressional  wave  portion  but  rather  on 
any  large  feature  of  the  seismic  wave  train  found  suitable  for 
computing  difference  in  time  of  arrival  to  field  emplaced  sensors. 

The  second  technique  relies  on  the  similarity  of  the  signals  from 
recoil  and  explosion  and,  therefore,  more  on  the  surface  wave  portions 
of  the  seismic  signature.  Both  techniques  for  location  require 
registration  or  pre-calibration  of  the  area  containing  the  energy 
mortar.  No  approach  was  found  which  could  utilize  seismic 
mortar  signature  data  singly  to  locate  position  of  enemy  fire,  i.e., 
additional  signal  arrival  time  information  is  needed  from  known 
source  locations  near  the  enemy  position  in  order  to  provide  a means 
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location.  It  remains  possible  that  these  sources  of  Information  could 
be  something  other  than  Impacts  if  locations  are  known.  To  obtain 
location  we  must  either  know  the  Individual  sensor  path  velocities 
to  the  enemy  gun  position  or  we  must  have  observed  Impulsive  signals 
from  a known  location  in  its  vicinity. 

3.1  HYPERBOLIC  LOCATION 

The  arrival  time  differences  of  seismic  surface  waves  to  possible 
pairs  of  seismometers  are  measured  and  locations  computed  by  the 
conventional  method  intersecting  hyperbolic  curves.  Assuming  a 
rectangular  coordinate  system  is  used  the  basic  equation  is  given  as 

Fj^j(x.y)  = [(x-x^)^  + (y-y^^)^]^'^^  - v^^/v^  [(x-x^)^  + (y-y^)^]^^^ 


- AT^j  v^ 


F^j(x.y) 


where  (x,y)  are  the  coordinates  of  the  seismic  source 

the  coordinates  of  the  two  sensors,  AT^^  is  the  arrival 

time  difference  between  the  two  sensors;  v, ,v.  are  the  two  source/ 

i J 

sensor  path  velocities.  The  only  difference  between  this  formulation 
and  the  conventional  one  used  for  sound  ranging  is  the  introduction 
of  individual  path  velocities.  Because  of  spatial  variations  in 
surface  geology  path  velocities  to  individual  spaced  sensors  cannot 
be  expected  to  be  the  same.  Having  estimates  of  the  velocities  these 
equations  may  be  solved  simultaneously  for  a sensor  triple. 

Path  velocity  estimates  are  obtained  by  timing  the  waves  generated 
by  return  fire  impacts.  These  path  velocity  estimates  are  then  used 
together  with  current  or  previously  measured  AT's  from  enemy  recoil  to 
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re-estimate  the  point  of  enemy  fire.  Return  fire  is  directed  to  this 
new  estimate.  Improved  estimates  of  path  velocities  are  obtained  by 
using  a Kalman  filter  [4]  to  provide  increasing  accuracy  to  the 
estimate  of  enemy  fire.  In  this  sense  the  location  algorithm 
adaptively  processes  the  sequential  data.  The  success  of  such  a 
processor  applied  to  acoustic  arrivals  would,  of  course,  be  dependent 
on  stability  of  atmospheric  conditions  which  effect  velocity. 

To  solve  these  equations  F(x,y)  we  can  linearize  them  by 
expanding  in  a Taylor  series  about  a trial  point  (*p»yp) 


F(x.y)  . F(yy^)  t f 


(x~x  ) + -r— 
P 9y 


P P 


(y-Vp) 

"P’^P 


which  is  a linear  form  of  F(x,y).  Having  estimates  of  the  path 
velocities  these  equations  may  be  solved  simultaneously  for  a sensor 
triple. 

An  alternative  method  of  solution  of  F^^(x,y)  = 0 and  F^j^(x,y)=0 
which  was  found  to  be  superior  is  obtained  by  formulating 


where  F . (x, y) =F  (x, y) =0  implies  the  surface  G(x,y)  is  at  a minimum, 
i J J 

Finding  such  a minimum  is  accomplished  by  the  method  of  steepest 
descent.  The  gradient  VG(x,y)  is  in  the  direction  of  maximum  change 
so  the  process  moves  by  trial  from  a starting  point  (x^,y^)  in 
the  direction  given  by  -VG(x,y).  After  a few  iterations  a good 
approximation  of  x,y  can  be  found  such  that 

Fij(x,y)  = F:ji^(x.y)  = 0 
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In  formulating  the  Kalman  filter  the  basic  state  equation 
consisted  of  the  path  velocities. 


Vu-V,  , + W,  , 
k k-1  k-1 


where  Vj^  is  the  state  vector  of  the  reciprocal  path  velocities  and 
W|^_^  is  the  state  Gaussian  noise.  Since  the  path  velocities  were 
assumed  to  be  constant  from  the  individual  sensors  to  the  locations  of 


enemy  fire  and  impacts  from  return  fire  Wj^  is  taken  as  zero  for  all 


iterations  k.  In  addition  to  the  state  equation  we  have  the  measure- 
ment equation: 

Z. 


Vk 


where  Z,  is  the  measurement  vector  of  observed  time  difference  of 
k 


arrivals  (AT's).  is  the  observation  matrix  containing  sensor 


impact  distances  d^j^  and  Is  the  Gaussian  measurement  noise  vector 


associated  with  Using  the  basic  time  difference  equation  for 

sensors  i and  j. 


AT 


ij 


'',1 


V . 
X 


The  matrix  is  2p-2  by  p with  form: 


“d2j^> 0 , 0 , . . . , 0 , 0 


0 » *^2k’  ~*^3k  > ® > ’ • ’ > ® ® 


0,...,0,...,0,d^p_^^j^,-dpj^ 


dik.-d2,0,0, . . . ,0,0 
0,  d'2j^,  -d'3j^,0,...,0,0 


0,...  ,0,...,0,d'^p_j^^j^,-d'pj^ 


where  the  d and  d*  are  distances  from  two  separate  shots.  In 
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applying  a Kalman  filter  data  sets  may  be  treated  sequentially  or 
simultaneously  with  the  same  resulting  estimate  of  Vj^.  Our  approach 
took  data  In  pairs  of  sets  to  avoid  the  difficulty  arising  from  the 
fact  that  there  are  p-1  AT's  and  p path  velocities. 

The  recursive  estimate  Vj^  of  the  state  vector  is  given  as: 


\ ■ \-i  * - \ Vi> 


where  is  the  Kalman  filter  matrix  and  minimizes  ^ I ^ 

The  filter  K^  is  obtained  from 

\ ‘ ‘■k-i  "k  <Vk-i  V * 


where  P,  = P,  , 
k k-1 


H,  P,  , and  R, 
k k-1  k 


covariance  of 


When  a new  measurement  Is  received  the  algorithm  has  already 

calculated  Pj^  and  ^ from  the  previous  iteration.  We  also  have 

by  calculation  of  the  sensor/impact  distances  using  the  current 

target  location  estimate.  P^  can  be  taken  as  the  identity  matrix 

and  V as  the  initial  estimate  of  the  path  velocities  in  order  to 
o 

start  the  recursive  process.  We  now  have  all  the  Information  needed 

to  obtain  a new  estimate  of  the  state  vector  V,  which  can  then  be 

k 

used  with  the  AT's  from  enemy  recoil  to  estimate  the  location 
coordinates  • 

Thus  this  algorithm  allows  one  to  obtain  an  updated  estimate 
of  the  location  of  enemy  fire  by  obtaining  successively  improved 
estimates  of  the  individual  seismic  path  velocities  to  that  point  of 
fire.  Since  the  mortar  recoil  signatures  are  highly  repeatable 
observation  of  successive  firing  can  only  help  improve  the  estimates 
of  AT  but  not  the  path  velocities.  Improvement  in  path  velocity 
estimates  are  obtained  by  timing  the  signals  from  known  Impact 
locations . 


’.'i 


3.2  METHOD  OF  RELATIVE  TIME  DIFFERENCE 

A second  algorithm  was  developed  based  upon  relative  comparison 
of  Impact  and  recoil  signals.  The  time  of  arrival  of  selsmlcs  from 
an  impact  are  the  same  as  a recoil  If  that  impact  occurs  at  the  same 
position  as  the  recoil.  Thus  the  observation  of  the  relative 
difference  In  arrival  times  between  an  Imoact  and  a recoil  can  be 
used  to  locate  one  event  with  respect  to  the  otner.  Furthermore,  If 
we  assume  that  the  paths  to  an  Impact  In  the  vicinity  of  a recoil  are 
unchanged  outside  of  path  length  and  the  velocity  of  the  seismic  waves 
are  constant  in  the  region  of  Impacts  and  recoil  then  location  of  the 
enemy  recoil  is  possible  if  the  approximate  position  of  the  impacts 
are  known.  Observation  of  the  difference  in  arrival  times  between  an 
impact  and  a recoil  is  accomplished  by  aligning  the  signals  at  one 
sensor  and  measuring  the  time  difference  between  the  others.  A 
minimum  of  four  sensors  are  required  to  obtain  an  estimated  location 
plus  the  wave  velocity  in  the  vicinity  of  the  known  impact  and  esti- 
mated recoil. 

The  solution  is  obtained  from  intersection  of  perpendiculars 
constructed  with  respect  to  lines  passing  through  each  sensor  and 
the  point  of  impact.  Figure  9 provides  a geometric  construction  of 
the  solution  for  four  sensors.  The  relative  AT's  as  calculated  above 
are  equal  to  the  actual  delays  (as  if  the  two  events  occurred  simul- 
taneously) plus  a constant.  Positions  on  respective  lines  are  now 
found  by  adding  in  an  arbitrary  constant  k and  moving  a distance 
(AT  + k)  • V from  the  point  of  impact.  At  each  position  a perpendi- 
cular P (or  circular  arc)  is  constructed. 

The  equation  of  the  family  of  lines  (perpendiculars  P)  parallel 
to  a line  which  passes  through  the  desired  location  and  which  is 
perpendicular  to  the  sensor  impact  line  is  given  by 

Y = M.x  + b . (k) 
j J 
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FIGURE  9.  GEOMETRIC  CONSTRUCTION  OF  SEISMIC  RANGING  ALGORITHM 
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where  is  the  Inverse  slope  of  a line  between  sensor  j and  the 
known  impact  (Xj.,Y^),  The  intercept  is  found  as 

b (k)  = (Y  - [AT,  + k]  V sin  6.)  - M, (X^  - [AT,  +k]  V cos  0,) 

J Ij  Jjlj  j 

where  0^  = tt/2  - TAN  The  intersection  of  such  perpendiculars 

as  determined  from  two  separate  sensors  form  a line  which  contains 
the  estimate.  The  intersection  for  sensors  i and  j has  coordinates: 


X*(k)  = (ATj  V sin  0j  - AT^  V sin  0^  + k V sin  0^ 


- k V sin  0.  - X^  - M,  AT,  V cos  6 . 

J i 3 I J j 3 

+ M,  AT,  V cos  0,  - k M,  V cos  0.  + k M,  V cos  0,)/(M.-M,) 

ii  i j 3 i i^'ji'^ 

Y*(k)  = Mj  X*(k)  + b^(k) 

The  line  of  intersections  is  given  parametrically  as  X*(k),  Y*(k). 

A second  line  of  intersections  generated  by  a second  set  of  two 
sensors  can  be  solved  simultaneously  with  the  first  to  produce  an 
estimate  of  location. 

If  arcs  are  generated  from  respective  sensor  positions  instead 
of  perpendiculars  a more  accurate  estimate  can  be  obtained.  This 
estimate  assumes,  however,  a uniform  known  velocity  throughout,  as 
was  the  example  case  of  Figure  9.  In  later  simulations  of  the 
method  of  arcs  certain  instabilities  resulted  which  did  not  occur 
when  using  perpendiculars.  For  this  reason  the  method  of  arcs  was 
discarded  in  favor  of  the  method  of  perpendiculars  in  conducting 
the  simulations. 
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While  there  are  similarities  between  this  approach  and  registra- 
tion as  used  for  acoustic  sound  ranging  there  are  important  differ- 
ences. First,  the  path  velocities  effects  are  limited  to  a region 
near  the  Impacts  and  recoils.  Secondly,  the  velocities  are  constant 
in  time  and  do  not  affect  the.  solution  by  an  great  degree. 

3.3  SIMULATION  OF  LOCATION  ALGORITHMS 

The  computational  algorithms  could  not  be  evaluated  in  forms 
of  seismic  field  data  recording  which  existing  prior  to  this  study 
program.  Simulation  was  chosen  as  a means  of  investigating  the 
location  capabilities  of  these  algorithms. 

Several  sources  of  error  are  present  in  computation  of  location 
estimates . 

1.  Sensor  location  error 

2.  Measurement  error  for  AT 

3.  Dispersion  and  aiming  error  for  return  fire 

4.  Errors  in  the  estimates  of  surface  wave  velocities. 

By  using  simulation  techniques  the  effect  of  each  of  the  errors 
can  be  examined  and  the  overall  performance  of  the  proceedures 
evaluated. 

Flow  diagrams  for  each  simulation  program  are  provided  in 
Figures  10  and  11.  Input  data  requirements  in  each  case  include  a 
rough  first  location  estimate  based  on  the  first  observed  enemy 
recoil.  In  actual  practice  a first  estimate  could  be  provided  by 
prior  calibration  of  the  area  or  use  of  air  wave  directionality. 

Error  sigmas  needed  for  simulation  were  essentially  unknown 
and  should  be  considered  carefully  when  examining  the  results.  Within 
each  simulation  program  the  error  sigma  were  used  to  generate  random 
values  for  each  parameter.  Representative  input  values  for  the 
simulation  programs  are  shown  in  Table  1.  Values  of  the  measurement 
errors  were  specified  at  100  m path  separation  distance.  For  paths 
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INPUT  DATA:  Sensor  Source  Geometry, 

Path  Velocities,  First  Velocity  and  Location, 
Estimates,  Error  Sigmas,  Number  of  Iterations 


Random  Values 
From  Error 
Distributions 


I Cycle  n Times^ 


H 

U 

< 

04 


o 

u 

u 

tfl 


Calculate  Apparent  Sensor  Locations 
Calculate  Recoil  Times  With  Error 


"1  Fire  Initial  Set  of  n Impacts  Around 
J Location  Estimate 


1 


Calculate  Impact  Times  With  Error 

Use  Kalman  Filter  to  Improve  Path  Velocities 


Use  Recoil  Times  With  Current  Location 
Approximation  to  Solve  Hyperbolic  Equations 
For  Each  Triplet  of  Sensors.  Average  All 
Solutions  For  New  Estimate 


I 


Random  Selection  of  Recoil  or  Impact 


IMPACT 


Return  Fire  to  Current  Location  Estimate 


CSTOP)- [ After  m Iterations  j 


FIGURE  10.  SIMULATION  PROGRAM  FOR  A HYPERBOLIC  FIXING 
LOCATION  ALGORITHM 
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If  No  Solution 
Adjust  K 


Calculate  Positions  Along  Each  Sensor 
Impact  Path  Which  Corresponds  to 
AT  + K, , AT  + K„ 

i Li 

Using  Sensors  S.  S.  Solve  For 

Intersection  of  Lines  or  Circles 

Determine  the  Line  of  Intersections 

Parametric  in  the  Sliding  Parameter 

K.  Repeat  for  Sensors  S.,  S,  Find 

r i+l 

the  Intersection  of  Two  Such  Lines. 


Average  the  n-2  Solutions  for  This  jth 

Iteration  F. 

3 

Combine  Averages  F^,  Fg,  ...  , F^  for 
the  Current  Approximation,  F Average 


After  m Iterations 


STOP 


FIGURE  11.  SIMULATION  PROGRAM  FOR  A SEISMIC  RANGING  LOCATION  ALGORITHM 
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TABLE  1 

SIMULATION  INPUT  PARAMETERS 


Parameter 

Values 

Comments 

Range  error  a 

K 

22.7m 

> 

81  mm  mortar  return  fire. 
Round  dispersion  errors. 

Deflection  error  o^ 

7.2  m 

Approximates  short  range 
artillery. 

Measurement  error 

5-40  ms 

Delay  timing  errors  at 
100  m. 

Sensor  Position  error  Og 

5-20  m 

Same  error  in  both  x and  y 
directions . 

Initial  Source  Location 
error 

100-1000  m ; 

Separation  distance  between 
the  true  location  and  first 
estimate. 

Path  Velocities 


AOO-1200  m/sec 


For  seismic  surface  waves 
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of  greater  separation  the  v;as  taken  to  be  proportionately  larger, 
i.e.,  0^(200  m)  »’  2o^(100  m)  . The  actual  relationship  between  o^ 
and  separation  is  unknown.  Fortran  computer  programs  for  each' 
simulation  algorithm  are  listed  in  Appendix  A.  These  programs  were 
built  up  on  an  IBM  370  system. 

The  effect  of  using  a Kalman  filter  is  demonstrated  in  Figure  12. 
Here  convergence  of  mean  path  velocity  errors  are  shown  as  a function 
of  the  number  of  impacts.  The  bars  indicate  the  range  of  path 
velocity  errors.  Beginning  with  a large  mean  path  velocity  error 
100  m/sec  (mean  path  velocity  = 600  m/sec)  and  a rather  large  measure- 
ment error  sigma  of  20  ms.  Eight  to  10  rounds  appear  sufficient 
for  adequate  convergence.  The  overshoot  which  appeared  at  4 rounds 
is  a characteristic  of  the  Kalman  filter. 

Typical  results  from  simulation  of  the  hyperbolic  algorithm  for 
various  values  of  Og  and  are  given  in  Table  2.  The  relationship 
between  the  location  estimate  error  and  the  average  velocity  error 
was  investigated  for  eight  successive  impacts  and  various  values  of 
O.J,.  With  increasing  scatter  the  points  plotted  in  Figure  13  show  that 
tile  location  error  increases  directly  with  the  average  velocity  sensor 
after  eight  impacts.  The  initial  average  path  velocity  error  was  for 
each  trial  in  excess  of  100  m/sec. 

Tne  input  sensor  geometry  consisted  of  five  to  seven  sensors  on 
a line  spaced  from  100  to  1000  meters.  The  spacing  used  in  a majority 
of  cases  was  500  meters.  The  position  of  enemy  recoil  was  located 
1500  to  2500  m from  the  center  of  the  baseline. 

Results  from  simulating  the  seismic  ranging  algorithm  have  been 
plotted  as  histograms  in  Figures  14  through  17.  Each  histogram 
describes  the  number  of  impacts  which  resulted  in  a location  estimate 
with  less  than  50  m error.  The  abscissa  shows  the  sequential  round 
number  with  percent  convergence  after  three  shots.  Each  histogram 
represents  200  independent  trials.  Individual  histograms  are  shown 
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Measurement  Error  = 2 
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Initial  Source  Error  707m 


12.  CONVERGENCE  OF  MEAN  PATH  VELOCITY  ERRORS 
USING  KALMAN  FILTER  SIMULATION  PROGRAM 


number  of  impacts 

FIGURE  14.  HISTOGRAMS  FROM  RANGING  LOCATION  SIMULATION.  INITIAL  LOCATION 
ERROR  EQUALS  200  m.  A BASELINE  OF  FIVE  SENSORS  WAS  USED. 


iNfilNG  LOCATION  SIMULATION.  Oj  = 10  ms 
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with  respective  error  assumptions.  In  Figure  14  histograms  resulting 
from  the  simulations  show  effect  of  increasing  error  sigmas  for  sensor 
location  and  time  measurement.  The  effect  of  increasing  initial 
location  error  can  be  seen  by  comparison  of  the  three  histograms  in 
Figure  15.  Figure  16  indicates  that  the  Initial  velocity  error  has 
little  effect  on  the  algorithms  location  capability.  Finally,  the 
effect  of  sensor  spacing  is  indicated  in  Figure  17. 

It  appears  from  these  simulation  results  that  the  error  in 
initial  location  is  most  sensitive  to  the  ultimate  ability  to  compute 
good  locations.  Sensor  spacing  and  initial  velocity  estimates  do  not 
greatly  influence  the  final  result.  The  effects  of  measurement  and 
sensor  location  errors  are  on  the  same  order  and  increasing  errors 
progressively  degrade  the  result.  In  order  to  achieve  acceptable 
results  from  the  location  programs  it  appears  that  measurement  sigmas 
should  not  exceed  10  ms  and  sensor  location  sigmas  should  not  exceed 
10  or  20  meters. 

These  results  are  from  simulation  and  as  such  cannot  provide  good 
evaluation  without  sufficient  field  data.  For  this  reason  the  study 
embarked  on  a field-measurements  program  in  order  to  provide  a more 
complete  basis  for  evaluation. 

3.4  DIRECTION  ANALYSIS  METHOD  FOR  ACOUSTIC  ARRIVALS 

Both  of  the  techniques  developed  require  a first  estimate  of 
location  to  start  the  process.  This  first  estimate  does  not  have 
to  be  highly  accurate  but  should  be  within  a few  hundred  meters  of  the 
target.  One  possible  method  which  could  be  used  to  obtain  this  type 
of  information  is  to  examine  the  directionality  of  the  seismic  and 
acoustic  signals  from  enemy  mortar  recoil  and  muzzle  blast.  Utilizing 
direction  (azimuth)  information  from  multiple  three  component  seismo- 
meter packages  a rough  location  can  be  calculated.  The  acoustic 
signals  potentially  have  greater  possibilities  for  source  direction 
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determinations  because  their  signals  are  much  larger.  The  following 
approach  based  on  signal  energy  was  developed  for  signal  direction 
determinations. 

The  direction  of  a wavefront  normal  is  the  direction  of  maximum 
energy  flow  in  the  horizontal  plane.  This  direction  can  be  found 
for  a bundle  of  energy  in  the  following  way.  Let  = L(t^)  and 
^i  ” digitized  signals  received  on  the  longitudinal 

and  transverse  seismometers,  respectively,  at  time  t^.  The  bundle  of 
energy  passing  during  the  time  from  t^  to  tj  will  be  proportional  to 


J 

E 

k-i 


on  the  longitudinal  component, 


E 


T 


j 

E 

k-i 


on  the  transverse  component,  and 


E ■ + total. 


If  the  seismometers  were  planted  at  an  angle  0 counter-clockwise  to  the 
way  they  were  planted  above,  the  signals  would  be 


\<0)' 


\^L^(0) 


cos  9 sin  6 


-sin  0 cos  9/  \ L 


and  the  energy  equations  would  be 


j 

E^(6)  I [T^(6)]2  and 
k“l 

E(0)  - E^(0)  + E^(0)  . 

The  direction  of  maximum  energy  flow  will  be  the  angle  0 from  the  longi- 
tudinal component  for  which  E^(0)/Ej^(0)  Is  a minimum.  In  other  words, 
we  rotate  the  seismometers  until  the  ratio  of  transverse  to  longitudinal 
energies  Is  a minimum.  Then  the  rotated  longitudinal  component  will  be 
parallel  to  the  Rayleigh  wavefront  normal. 

A closed  form  of  the  angle  0 for  which  E^(0)/Ej^(0)  is  minimal  can 
be  derived  as  follows.  We  first  note  that  E(0)  » E,  i.e. , the  total 
energy  in  the  horizontal  plane  is  not  dependent  on  the  way  the  seismometer 
Is  planted: 


E(6)  = E^(0)  + E^(0) 


",  ° E [T,(9)]2  + E [L  (0)]2 

k-i  k-i 
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j 

- E m + L2)2 


j j 

" r T2  + E L2 
k»i  k-1  ^ 


- E 

We  now  can  see  that  E^(e)/Ej^(e)  is  minimal  when  E^^Ce)  is  a maximum: 

E^(e)  + Ej^(e)  - E - const, 
or 

Et(0)/Ej^(6)  ■ const/Ej^(6)  - 1 . 

The  right-hand  side  is  clearly  minimal  when  Ej^(0)  is  maximal.  To  find  the 
desired  angle  we  find  0 for  which  E^(0)  is  a maximum. 
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_d 

d6 


J 2 

z [L.(e)r 

k-i 


J 

- 2 s L (e) 

k-i 


J_ 

de 


[Lj^(e)l 


Setting  this  equal  to  zero. 


0 . z v«>  i 

k-i 


- E (-sine  + COS0  Lj^')  (-sine  Tj^  + cos6  Lj^) 
k-i 


- I (-sine  Tj^  + cose  Lj^)  (-cose  Tj^  -sine  Lj^) 
k-i 


- sine  cose  Z T?  + (sin^e  - cos^-0)  Z I.  L -sin9  cose  Z 
k-i  ^ k-i  k-i 


J , ^ ^ 2 . 

Let  A - Z t2  B - Z T L , C - Z and 

k-i  k-i  ^ k-i 


divide  by  -sine  cos0  : 


- A - 


sin^e  - cos^e 
sine  cose 


0 


This  can  be  written  as 


tane 


1 

tane 


9 
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or  as 


tan^e  + — -■  tan6  - 1 - 0 , 


Of  the  two  solutions,  the  0 which  makes  E^O)  larger  Is  the  angle  between 
the  longitudinal  component  and  the  Rayleigh  wavefront  normal. 


A 

SYSTEM  DESIGN  CONSIDERATIONS 

The  system  will  be  based  upon  a distributed  field  of  seismometers 
which  partially  surround  the  enemy  weapons.  Given  the  arrival  times 
of  several  mortar  recoil  signals  the  problem  reduces  to  computation  of 
locations  from  sets  of  three  or  four  sensors  by  the  method  of  inter- 
section of  hyperbolic  curves.  Discrimination  of  the  seismic  recoil 
signals  from  shell  impact  signal  is  necessary  when  multiple  events 
can  occur  with  rapid  succession.  In  this  case,  arrival  must  be  sorted 
out  according  to  source.  The  shell  impacts  from  friendly  return  fire 
can  be  used  to  provide  calibrating  path  velocities  and  can  be 
incorporated  into  a "seismic  on  seismic"  location  estimate.  At 
present,  the  seismic  mortar  system  is  envisioned  to  consist  of; 

(1)  several,  three  orthogonal  component  sensor  units  equipped  with 
signal  processing  electronics  and  data  transmission  capabilities, 
and  (2)  a central  (base  operated)  receiver,  processor,  location 
computer,  and  display  unit.  A block  diagram  is  shown  in  Figure  18. 

The  first  three  blocks  constitute  the  remote  battery-powered 
sensor  package.  The  necessary  data  from  this  sensor  could  be 
digitized  and  transmitted  by  radio  frequency  (RF)  link  to  a small 
processing  computer  controlled  by  the  system  operator. 

In  the  present  concept  the  sensor  package  should  be  small,  light 
weight,  low  cost,  and  perhaps  expendable.  The  unit's  principal  func- 
tions will  detect,  identify,  and  report  impacts  and  recoils.  The 
unit  will  likely  not  be  able  to  discriminate  weapon  type  or  source. 

The  sensor  will  become  activated  by  a basic  amplitude/time  threshold. 
With  threshold  activation  incoming  analog  signals  will  be  digitized 
and  stored  for  data  transmission  and  identification.  Perhaps  the 
best  available  means  for  identification  can  be  made  by  observing  the 
air  wave  arrival.  The  logic  would  say  an  event  has  occurred  within  a 
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FIGURE  18.  SEISMIC  MORTAR  LOCATION  SYSTEM 
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proper  time  lag  of  the  threshold  seismlcs  and  which  has  suitable 
amplltude/f requency/duratlon  characteristics.  Thus  the  seismic 
activation  is  confirmed  as  that  from  an  impact  or  recoil  by  presence 
of  an  air  wave. 

Data  burst  transmission  should  begin  as  soon  as  possible  after 
the  air  wave  arrival.  Depending  on  RF  data  transmission  capabilities 
signals  may  have  to  be  staggered  to  prevent  overlapping  transmissions 
from  several  activated  sensors. 

Without  special  provisions  the  minimum  time  spacing  between 
events  may  be  as  much  as  five  seconds.  If  a push  down  stack  or  some 
similar  technique  is  employed  this  time  could  be  reduced  to  as  little 
as  two  or  three  seconds.  The  lower  limit  would  be  predicated  by 
desired  seismic  trace  length. 

When  rapid  fire  occurs  seismic  signatures  will  overlap  in  time 
which  may  destroy  their  usefulness  for  timing  information. 

The  remote  battery-powered  sensor  must  be  able  to  generate  an 
electric  signal  corresponding  to  ground  velocity  of  adequate  signal- 
to-noise  ratio  (S/N)  in  the  range  of  seismic  frequencies  (10-40  Hz) 
that  contain  mortar  location  information.  A small,  light-weight, 
compact  seismometer  such  as  the  Sparton  device,  with  characteristics 
shown  in  Figure  19  would  provide  sufficient  sensitivity.  A preampli- 
fier with  a gain  of  about  4000,  not  a difficult  requirement,  would 
then  be  necessary  to  increase  this  maximum  seismometer  output  to  1 
volt,  an  amount  appropriate  for  the  succeeding  circuitry.  A 
seismometer  with  greater  sensitivity  could  be  used  but  it  would 
probably  be  larger,  heavier,  and  more  fragile;  so  it  seems  more 
reasonable  to  use  electronic  amplification  instead.  At  the  above 
stated  signal  levels  the  S/N  ratio  would  not  be  degraded  by  electronic 
noise.  The  amplifier  could  be  a metal  oxide  semiconductor  (MOS) 
integrated  circuit  (IC)  that  could  operate  from  small  batteries  for  90 
days  without  difficulty.  Total  system  power  requirements  will  be 
discussed  later  in  this  section. 
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Because  there  are  many  frequency  components  generated  in  a 
battlefield  situation,  the  low  pass  (LP)  and  high  pass  (HP)  filters, 
are  necessary  to  remove  unnecessary  and  detrimental  frequencies. 

In  this  manner,  the  S/N  ratio  of  the  electrical  signal  may  be 
enhanced.  Adjustment  of  amplifier  gain  and  the  placement  of  LP  and 
HP  cutoff  frequencies  will  be  possible  to  adapt  the  sensor  package 
to  that  geology. 

The  remaining  functions  in  the  sensor  are  the  threshold  circuit, 
A/D  conversion,  and  an  RF  data  link.  It  is  not  an  effective  use  of 
battery  power  to  transmit  data  continually  since  useful  data  only 
occurs  during  and  immediately  after  each  shot.  Therefore,  a circuit 
is  necessary  to  determine  when  the  seismic  signal  has  exceeded  a pre- 
set threshold  and  then  the  transmitter  is  turned  on  for  a short  time 
to  relay  the  information.  This  approach  is  also  desirable  from  the 
security  standpoint  because  short  bursts  of  RF  are  more  difficult 
to  monitor  and  detect  by  the  enemy.  Additionally,  RF  bursts  will 
occur  during  mortar  acoustic  reports  which  means  that  acoustic  indi- 
cations from  the  enemy's  monitor  receiver  will  be  masked  by  the 
external  muzzle  and  impact  sounds. 


The  central  processing  unit  contains  the  equipment  to  receive 
and  detect  the  RF  signals.  The  source  of  the  digital  data  transmission 
can  be  identified  by  a sensor  code.  Digital  samples  will  be  stored 
in  memory  locations  known  to  the  computer  program.  Storing  of  all 
sensor  outputs  will  begin  when  the  first  RF  transmission  occurs. 

To  store  only  useful  data,  the  digitizing  will  continue  for  that  shot 
or  series  of  shots  only  sufficiently  long  enough  to  collect  the 
seismic  wave  data. 

Only  a small  fraction  of  program  time  will  be  used  to  acquire  new 
data.  The  remaining  block  of  Figure  18  represents  functions  performed 
by  the  processor  during  most  of  the  program  time.  Actually,  the 
main  program  will  be  interrupted  relatively  infrequently  to  gather 
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new  data.  The  keyboard  and  display  may  also  be  serviced  on  an 
Interrupted  basis  auch  as,  when  a command  changes  or  when  the  visual 
display  needs  updating. 

4.1  POWER  CONSUMPTION 

Concerning  the  sensor  package  only,  the  power  requirements  of  the 
RF  transmitter  are  the  moat  difficult  to  establish  and  may  be  the  most 
severe  if  high  transmitter  power  Is  needed.  For  this  range  and 
frequency  over  water,  one-half  watt  of  RF  power  input  to  the  transmitter 
should  be  sufficient,  assuming  the  usual  military  band  of  buoy 
frequencies  of  about  170  megahertz  is  used.  To  be  ultra-conservative 
and  for  land-based  use  if  one  assumes  that  a transmitter  of  5 watts 
input  power  is  required  and  that  a single-event  transmission  takes 
2 seconds;  then  5 time  2 or  10-watt  seconds  are  needed  for  each 
event.  The  energy  density  of  an  advanced  nickel-zinc  battery  is  about 
40  watt-hours  per  pound  or  144,000  watt-seconds  per  pound.  This  means 
that  a one-pound  battery  should  transmit  for  14,400  events.  If  a 
90-day  life  is  necessary,  then  a firing  rate  of  90  x 24  x 60  or  33,600 
minutes  per  month  divided  by  14,400  possible  transmissions  or  about 
2 shots  per  minute  can  be  tolerated.  This  is  the  rate  for  a group 
of  weapons. 

The  remaining  electronics  in  the  sensor  package  have  minimal 
power  requirements.  assuming  10  operational  amplifiers  are  needed, 
with  each  using  1000  ywatts,  a conservative  estimate  for  low-power 
devices,  10,000  uwatts  times  2500  hours  in  90  days  is  25  watt  hours, 
less  than  a one-pound  battery.  This  amount  of  power  must  be  adequate 
with  the  sensor  either  operating  as  an  RF  device  since  these  circuits 
are  operating  continuously.  It  should  be  understood  that  this  estimate 
of  required  battery  power  is  very  conservative.  A 0.5  lb  battery  will 
buy  10  days  of  operation  and  3 km  transmission  range  which  is  a more 
practical  estimate. 


t'ORMtllLV  WILLOW  RUN  I.ARORATORIKR.  THI  UNIVIRRITV  OR  MICHIOAN 


All  sensors  will  send  their  data,  via  RF  to  a central  location. 

If  this  digital  processlr  Is  of  MOS  design,  power  requirements  will  be 
low  especially  If  an  electro-mechanical  disc/drum  memory  Is  not 
necessary.  The  addition  of  the  LED/llquld  crystal  display  and  control 
functions  completes  the  processor  with  the  total  power  requirements 
estimated  at  10  watts  without  any  mechanical  components.  This  10 
watts  could  be  provided  by  a small  5-lb  rechargeable  battery  unit  good 
for  one  or  two  days  of  operation.  If  the  batteries  are  of  the  quick 
recharge  type  such  as  lead-cobalt  they  could  be  used  until  recharging 
was  necessary  whep  a vehicle  engine-mounted  generator  could  quickly 
recharge  them  for  further  operation.  The  volume  and  weight  of  the 
sensor  units  and  central  computer  are  such  that  two  or  three  men  could 
carry  the  system.  The  total  weight  of  the  system,  including  10  sensor 
packages,  would  be  on  the  order  of  75  lbs.  This  estimate  assumes  a 
10-day  sensor  battery  life.  A total  system  cost  (on  a production 
basis)  is  estimated  to  be  $5K  to  $6K.  An  estimate  of  component 
sizes,  weights,  and  costs,  based  on  prior  experience  with  electronics, 
are  given  in  Table  4. 

In  summary,  the  system  just  described  should  not  be  difficult  to 
design  and  deploy.  There  are  no  components  approaching  state-of-the- 
art  limitations.  There  have  been  many  types  of  buoys  built  previously 
and  quite  some  time  ago.  In  the  digital  processing,  some  custom  LSI 
design  may  be  necessary  but  that  too  has  been  done  before,  though  more 
recently.  The  buoys  should  not  contain  many  components  so  that  failure 
rate  during  battery  life  could  be  made  low  by  selection  of  components 
and  adequate  specification  and  testing.  Plug-in  spare  units  or  circuit 
cards  could  be  used  to  increase  the  MT3F  (mean  time  between  failures) 
of  the  central  processing  gear  since  there  are  many  more  components 
involved  than  in  the  buoys.  The  most  fragile  element  of  the  entire 
system  would  probably  be  the  magnetic  drum  or  disc  to  head  interface 
if  that  type  of  storage  were  used. 


TABLE  3 


Eat|.inated  Component  Volume,  Uelghc,  and  Cost  of  the 
Seismic  Mortar  Location  System 


Component 

Size 

Weight 

Est.  Cost* 

(cu.  In.) 

, (lb) 

($) 

Seismometer 

1 

0.5 

5 

Amplifier  and 

8 

0.5 

20 

Filters 

Threshold,  Switch 

1 

0.5 

5 

and  Line  Driver 

A/D  Conversion 

30 

0.5 

50 

RF  Transmitter 

4 

0.5 

20 

Battery 

6 

0.5 

(10  day) 

10 

Enclosure 

64 

1.0 

10 

Total 

4x4x4 

4.0  lbs/ 
unit 

70/unit  * 

Data  Processing 

Receivers 

0.12 

3.0 

1000 

Central  Processing 

0.25 

10.0 

2000 

Unit 

Digital  Storage 

0.50 

2.0 

500 

Display,  Control,  0.50 

Batteries,  and  Enclosure 

20.0 

1000 

Total 

1.37 

13  X 13  X 13 

35  lbs/ 
unit 

4500/unit 

*Based  on  100  units 
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FIELD  MEASUREMENTS  PROGRAMS 


In  order  to  evaluate  the  use  of  seismic  data  for  mortar  location 
a series  of  field  measurements  was  initiated  during  the  second  phase 
of  the  project.  Seismic  mortar  and  explosion  data  were  collected  at 
two  sites.  Twenty  Nine  Palms  and  Camp  Pendleton,  California  during 
February  and  March  of  1975. 

The  object  in  this  field  exercise  was  to  simulate  as  near  as 
practical  those  battlefield  conditions  of  mortar  ranging.  To  do 
this  we  needed  to  make  recordings  of  enemy  recoil  and  friendly  shell 
bursts.  Enemy  mortar  fire  can  be  represented  by  firing  conventional 
81  mm  mortars.  It  was  not  practical  to  observe  the  effect  of  friendly 
fire  directly.  To  impact  the  vicinity  of  a mortar  firing  position 
would  have  required  setting  up  in  an  impact  zone  which  is  generally 
not  permitted.  Live  impacts  have  the  additional  problem  that  each 
must  be  surveyed  to  an  accurate  position. 

During  a previous  program  which  included  field  work  at  Ft.  Sill, 
Oklahoma  it  was  determined  that  105  mm  shell  impacts  could  be 
simulated  seismically  with  a TNT  detonation.  This  technique  was 
used  exclusively  throughout  the  field  work  in  order  to  record  the 
effect  of  friendly  fire. 

Reflecting  the  requirements  of  the  derived  algorithms  specific 
objectives  of  the  field  program  included  collection  of  data  which  could 
evaluate  the  similarity  of  impacts  and  recoils  as  a function  of 
separation  distance.  In  addition  we  had  hoped  to  evaluate  the 
concept  of  firing  into  an  area  (simulated  by  explosions)  for  purposes 
of  calibration  of  path  velocity. 
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5.1  FIELD  MEASUREMENT  PROGRAM  - TWENTY  NINE  PALMS  MARINE  BASE 

The  field  measurements  at  Twenty  Nine  Palms  were  made  at 
locations  near  Deadman's  Lake  and  Camp  Wilson  (see  map  in  Figure  20). 
The  topography  was  flat  in  this  area.  Surface  material  in  this  area 
consisted  of  unconsolidated  rock  and  sand  with  an  underlying  bedrock 
which  consisted  mostly  of  granite.  Depths  to  the  bedrock  material 
was  guessed  to  be  a few  hundred  feet.  After  preliminary  examination 
of  seismic  recording  was  made  in  this  area,  it  was  felt  that  most  of 
the  seismic  energy  travelled  along  paths  in  the  upper  unconsolidated 
material . 

A series  of  nine  three-component  seismometer  packages  were 
instrumented  along  a line  parallel  to  and  located  300  meters  inside 
the  southern  boundary  of  the  BRAVO-1  impact  area.  Mortar  firing 
positions  were  located  along  a 2.5  kilometer  line  of  0.5  kilometer 
intervals  and  perpendicular  to  the  sensor  base  line.  The  main  firing 
line  ran  from  the  center  of  the  sensor  line  (5th  sensor  position)  to 
a position  near  Camp  Wilson.  Two  pound  TNT  shots  were  to  be  fired 
in  a cluster  pattern  around  each  mortar  firing  position  simulating 
friendly  impacts.  A plan  of  this  field  geometry  is  shown  in 
Figure  21.  A standard  81  mm  mortar  was  used  to  fire  rounds  over  the 
sensor  line  into  the  BRAVO-1  impact  area.  All  explosion,  seismometer, 
and  firing  positions  were  surveyed  by  the  base  surveying  team. 

5.1.2  INSTRUMENTATION 

Figure  22  depicts  the  instrumentation  which  was  used 
in  digitally  recording  the  seismic  data  in  the  field.  This  system, 
housed  in  a van,  is  capable  of  recording  up  to  15,000  samples  per 
second  with  a dynamic  range  of  84dB,  14  bits  per  sample,  full 
frequency  resolution  being  limited  to  75  Hz  by  the  antialiasing 
filters.  The  controlling  computer  was  an  Interdata  Model  7/16 
processor  with  16,384  eight  bit  bytes  of  1,000  ns  core  memory. 
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Existing  peripheral  equipment  Included  an  ERIM  built  1.2  megacycle 
clock,  an  Analogic  aeries  AN5800  A/D  and  D/A  converter,  an  ASR 
model  33  teletypewriter,  an  Interdata  Hexadecimal  I/O  panel  and  a 
Kennedy  9000  digital  tape  recorder. 

Nine  three-component- othogonally  mounted  Hall  Sears  model  HS-10 
seismometers  were  used.  These  velocity  transducers  have  a resonant 
frequency  of  2 Hz  and  a damping  factor  of  approximately  0.64.  They 
have  been  calibrated  on  ERIM's  shaketable.  Each  Instrument  has  a 
response  to  motion  perpendicular  to  its  axis  which  Is  30  dB  down  from 
Its  response  to  on-axls  motion.  Figure  23  Is  a typical  amplitude 
response  curve  for  these  Instruments.  Each  seismometer  was  connected 
to  the  field  recording  van  using  standard  WD-1  communication  wire. 

Each  signal  was  amplified  at  the  seismometer  to  reduce  the 
percentage  of  induced  noise  in  the  line  carrying  the  signal  to  the 
recording  van.  Ithaco  model  9121-73  seismic  amplifiers,  powered  by 
automotive  type  12  v.  lead-acid  batteries,  were  used  here.  They 
can  be  adjusted  from  -12  to  96  dB  In  6 dB  intervals  and  are  capable 
of  providing  an  output  signal  of  5 v.  rms.  Their  low  and  high  cutoff 
frequencies  are  0.3  Hz  and  1,000  Hz,  respectively,  and  they  generate 
2 pv  equivalent  input  noise.  Zero  times  were  generated  for  each 
event  by  radio  transmission  of  the  immediate  seismic  impulse. 

Research  at  this  laboratory  has  shown  that  the  signals  of  interest 
for  long  range  detection  of  recoil  have  a spectral  peak  near  15  Hz 
with  energy  possibly  as  high  as  20  Hz.  Therefore,  filters  with  a 
response  as  shown  in  Figure  24  were  placed  in  the  system  before  the 
data  was  digitized  to  prevent  the  aliasing  of  higher  frequency  energy 
into  the  wanted  data  band.  Since  the  sampling  frequency  was  300  Hz, 
input  data  to  the  system  had  to  be  attenuated  an  amount  equal  to  the 
dynamic  range  of  the  system  at  the  Nyquist  frequency  or  150  Hz.  Tlijs 
condition  was  Insured  by  the  combination  of  the  use  of  the  filters 
and  the  character  of  the  mortar  recoil  seismic  signals. 
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FIGURE  23.  TYPICAL  RESPONSE  OF  GEOSPACE 
MODEL  HS-10  AS  CALIBRATED  ON 
ERIM'S  SHAKETABLE;  DAMPING 
FACTOR  OF  APPROXIMATELY  0.64 
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FIGURE  24.  ANTI- ALIAS  FILTER  RESPONSE.  Output 
voltage  is  dB  relative  to  OdB  input 
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The  filter  outputs  were  sequentially  sampled  and  digitized  with 
a skew  of  1/85,000  second  between  adjacent  channels.  The  digitizing 
rate  was  300  samples  per  second.  Full  scale  A/D  conversion  is  -5v  to 
+5v.  Amplifier  gains  were  set  so  that  the  input  seismic  signal  plus 
noise  were  within  6 dB  of  full  scale  without  over-recording.  The 
samples  were  stored  in  a buffer,  under  software  control,  and  then 
written  on  digital  tape. 

The  quality  of  the  data  was  monitored  continually.  Oscilloscopes 
were  used  at  the  input  panel  during  an  event  and  hard  copy  traces 
were  produced  after  each  event  on  a Seimens  6 channel  oscillograph. 

5.1.3  DATA  COLLECTION 

The  seismic  data  collected  at  the  Twenty-Nine  Palms 
site  was  generally  quite  good  but  not  totally  free  of  problems.  Bad 
weather  and  uncontrollable  delays  prevented  the  field  team  from 
collecting  all  data  as  originally  planned.  Mortar  firings  were 
completed  at  the  0.5,  1.0,  1.5,  and  2.0  km  sites  using  standard 
HE  rounds.  Firings  were  made  at  charge  7 or  8 and  at  elevations  to 
yield  a 3.0  or  3.5  km  range.  All  mortar  firings  were  made  by 
Camp  Pendleton  infantry  mortar  crews.  Because  of  delays  and  high 
wind  conditions  no  data  was  collected  at  the  2.5  km  site  and  mortar 
signatures  from  the  2.0  and  0.5  site  were  noisey  and  not  suitable 
for  further  analysis.  Time  limitations  prevented  the  crew  from 
making  all  the  explosion  shots  as  planned.  However,  complete  sets 
were  obtained  for  the  1.0  and  1.5  km  sites.  A total  of  40  explosions 
and  38  mortar  firings  were  recorded  at  the  Twenty  Nine  Palms  site. 

The  stars  next  to  the  surveyed  locations  of  the  field  geometry  map 
of  Figure  22  show  positions  for  which  explosions  were  detonated. 

All  explosions  were  fired  by  base  EOD  personnel. 
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5.1.4  DATA  QUALITY 

Data  quality  was  very  good  for  a major  portion  of  the 
events  recorded  except  during  windy  conditions.  Signal-to-noise 
•levels  were  generally  high  due  to  the  quiet  nature  of  this  area. 

Some  Instrumentation  problems  did  occur  for  a couple  of  channels  at 
sensor  one  and  sensor  nine  later  In  the  field  session.  The  problem 
was  not  detected  since  other  sensor  packages  were  being  monitored. 

5.2  FIELD  MEASUREMENT  PROGRAM  - CAMP  PENDLETON,  CALIFORNIA 

A second  group  of  seismic  recordings  were  made  at  Camp  Pendleton 
firing  ranges.  This  area  has  a more  rugged  terrain  and  high  seismic 
noise  levels  generated  by  base  activities  which  presented  a sharp 
contrast  to  conditions  found  at  Twenty  Nine  Palms. 

It  was  more  difficult  to  set  up  in  this  area  because  of  the 
hilly  terrain  and  the  constraints  imposed  by  the  firing  range  layout. 

A map  of  the  Camp  Pendleton  site  area  is  shown  in  Figure  25  with  the 
placement  locations  of  the  nine  seismometer  packages  in  the  X-ray 
impact  area.  Mortar  firings  were  made  from  ranges  108,  108A,  and 
111.  Range  108A  was  located  approximately  halfway  between  range  108 
and  range  115.  Figure  26  shows  the  approximate  location  of  the 
explosion  shots,  which  simulate  impacts  in  the  area  of  each  mortar 
firing  position.  Because  of  the  terrain  we  were  unable  to  lay  them  out 
in  a uniform  pattern. 

The  area  surface  materia]  was  felt  to  consist  largely  of  hard- 
rock  with  a shallow  layer  of  surface  soil  material.  On  steep  slopes 
rock  outcropping  was  frequent.  The  recording  area  exhibited  several 
hundred  feet  of  relief.  Basically  there  was  a single  ridge  system 
separating  the  seismometer  positions  located  in  an  outwash  from  the 
firing  positions  on  range  108.  From  range  111  this  separation 
amounted  to  a two  ridge  system. 
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All  Instrumentation  was  the  same  as  that  described  earlier. 

5.2.1  DATA  COLLECTION 

A total  of  53  two  pound  shots  and  27  mortar  rounds  were 
fired  at  the  Camp  Pendleton  site.  Many  of  these  events  were  dupli- 
cate measurements  which  were  found  necessary  because  of  the  high 
general  noise  conditions  and  frequent  Interference  by  tanks,  trucks, 
artillery,  and  helicopters.  Noise  levels  at  Camp  Pendleton  were 
approximately  20  dB  higher  than  those  which  were  common  at  the  Twenty 
Nine  Palms  site.  Mortar  firings  were  made  at  the  standard  positions 
at  each  range.  Rounds  were  fired  at  charge  7 and  at  high  elevations. 
Explosions  were  detonated  In  the  vicinity  of  each  mortar  firing  pos- 
ition. Locations  surveyed  by  a Marine  Corps  support  team  are  shown 
In  Figure  26  with  starred  position  Indicating  those  which  were 
utilized  during  the  recording  program. 

5.2.2  DATA  QUALITY 

Except  In  one  or  possibly  two  Instances  seismic  signals 
from  81  mm  mortar  recoil  could  not  be  observed  at  any  sensor.  The 
exceptions  are  of  poor  quality  and  uncertain  because  of  the  large 
number  of  inter f erring  sources.  The  two  pound  explosions  were 
generally  observed  selsmically  out  to  the  fifth  or  sixth  sensor 
position  but  lost  in  the  noise  at  greater  ranges.  Many  of  the 
recordings  were  rendered  useless  because  of  interferring  events. 
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DATA  ANALYSIS  AND  RESULTS 

Our  approach  In  analyzing  the  field  data  was  to  first  examine  the 
data  visually  and  eliminate  those  events  which  were  not  considered  use- 
ful because  of  their  high  noise  levels  on  instrumentation  problems. 

Basic  processing  Including  Identification  events  zero  times » and  corre- 
lation parameters  was  completed  on  a PDP-8  special  purpose  computer 
system.  Seismic  and  acoustic  travel  times  were  calculated  using  first 
arrival  and  peak  correlation  sample  locations.  Actual  travel  times 
or  unnormalized  relative  times  for  each  event  along  with  shot  and  sensor 
position  information  served  as  Input  to  the  location  computer  code. 

These  data  were  not  computer  analyzed  for  purposes  of  discrimination  of 
recoil  from  explosion  events. 

6.1  RESULTS  OF  VISUAL  EXAMINATION 

The  field  measurement  program  produced  a total  of  4,266  records 
each  with  a seismic  and  acoustic  signature.  It  was  considered 
impractical  to  process  every  possible  record.  Visual  selection  of  a 
subset  of  data  based  upon  the  following  criteria. 

1.  Where  events  were  repeated  the  event  with  best  signal-to-noise 
was  selected. 

2.  The  largest  seismic  and  acoustic  signals  were  observed  on  the 
vertical  channel.  Also  greatest  correlation  between  shots 
was  observed  on  the  vertical  oriented  seismometer.  For  these 
reasons  only  the  vertical  channel  was  processed  with  the 
exception  of  directional  analysis. 

3.  Since  the  signatures  from  mortar  recoil  were  repeatable 
only  two  or  three  events  from  each  firing  site  were  selected 
for  processing. 

The  events  analyzed  include  most  of  those  which  have  circles  in 
Figures  21  and  26.  A list  of  the  survey  coordinates  of  all  these 
events  is  given  in  Appendix  B . 
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The  frequency  content  of  the  seismic  signatures  was  examined 
using  variable  Kronkite  analog  filters  in  conjunction  with  the 
digital  recording  system.  Figure  27  shows  example  passbands  of  an 
explosion.  Note  that  only  slight  noise  reduction  is  obtained  after 
using  the  filter  passband.  The  principal  frequency  components  of  the 
seismics  occur  near  twenty  hertz  with  most  of  the  energy  at 
frequencies  less  than  forty  hertz.  Unfortunately,  some  of  the 
principal  spectral  components  of  wind,  vehicles,  and  aircraft,  l.e., 
in  this  frequency  band  (10-40  Hz).  It  was  concluded  that  signal-to- 
noise  improvement  could  not  be  greatly  enhanced  by  simple  passband 
techniques.  Visual  examination  of  the  three  component  data  indicated 
two  features: 

1.  the  surface  waves  on  the  vertical  channel  were  usually  the 
largest  seismic  signal;  and 

2.  the  surface  waves  showed  a great  deal  of  similarity  from  shot 
to  shot.  Figure  28  shows  examples  of  the  three  component  analog 
records  from  two  explosions  which  were  fired  at  different  positions. 
Figure  29  shows  the  vertical  axis  seismics  from  two  explosions 
separated  by  125  meters  as  received  at  each  of  the  seismometer  array 
stations.  The  events  overlay  one  another  to  demonstrate  the  similarity 
of  signals.  Also  note  that  these  signals  appear  more  similar  between 
events  than  as  between  stations. 

6,2  TRAVEL  TIME  CALCULATIONS 

Acoustic  and  seismic  travel  times  for  the  first  arriving 
signals  were  calculated  by  first  locating  the  zero  time  and  picking 
the  first  arrival  to  the  nearest  digital  sample.  For  acoustic  signal 
the  resulting  travel  time  is  similar  to  that  which  would  be  obtained 
by  thresholding.  Due  to  the  variances  in  seismic  record  and  noise 
levels  simple  thresholding  would  prove  very  inaccurate.  Timing  was 
made  on  the  basis  of  common  feature  recognition  along  the  linear 


69 


>.i  iiilllliliM  W'jlfi^iiMIP 


70 


aUMBMIB^gg 


rORMIRLY  VYILLOW  RUN  LA»OR*TORI»S,  THE  UNIVERBCTY  Ol'  MICKIOAN 


2ip 


FO»MI(»».V  WILUOW  «UN  L*»OHATO«IE».  THE  UNIVEBEITV  OF  MICHIOAN 


array.  Using  seismic  first  arrivals  presented  some  difficulty 
because  this  signal  is  usually  not  the  largest  and  is  often  obscured 
by  background  noise.  Some  travel  time  Information  is  contained 
in  Appendix  B.  Crops  correlation  analysis  was  also  used  to  gain 
time  of  arrival  and  is  discussed  in  a later  section.  Air  coupled 
Rayleigh  waves  were  observed  as  distinctive  pre-acoustic  wave 
arrivals.  They  were,  however,  not  consistently  observed  and  therefore 
not  considered  for  analysis. 

6.3  SOURCE  DIRECTION  DETERMINATIONS 

As  previously  discussed  each  location  algorithm  requires  an 
initial  rough  location  estimate  based  on  the  first  observations  of 
enemy  fire.  One  possibility  is  to  obtain  signal  directivity  from 
each  three  component  sensor  by  measuring  the  relative  signal  (mostly 
acoustic)  energy  on  each  horizontal  channel.  A computer  program  was 
written  based  upon  the  computational  method  discussed  in  section 
A total  of  15  events  from  the  Twenty  Nine  Palms  data  were  used  to 
test  this  approach  which  had  appeared  encouraging  from  a few  simple 
in-house  laboratory  experiments.  Approximate  angular  deflections 
were  estimated  to  the  nearest  five  degrees.  Figure  30  shows  the 
typical  output  of  the  processor  with  the  top  trace  representing  the 
estimated  angular  deflection  with  time  and  the  bottom  trace  showing 
the  energy  rate.  In  order  to  obtain  an  estimate  the  top  trace  is 
weighted  by  the  bottom  trace.  The  results  are  summarized  in  Table  5 
They  were  not  very  encouraging.  Overall  the  method  could  pick  90 
degrees  window  which  contained  the  correct  azimuth  78  percent  of  the 
time.  One  possible  source  for  the  error  could  be  corregations  in  the 
acoustic  wavefront.  While  this  method  did  not  prove  successful  further 
work  is  needed  to  fully  investigate  this  possibility.  Other 
possibilities  of  obtaining  initial  location  estimates  Include,  of 
course,  using  acoustic  and/or  seismic  data  with  velocity  guesses 
and  the  hyperbolic  method. 
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FIGURE  30.  PLOT  OF  THE  TIME  VARYING  ANGULAR  DEFLECTIONS  AS 
DETERMINED  BY  THE  COMPUTER  ALGORITHM  FROM  THE 
ACOUSTIC  SIGNATURE.  SHOT  307.  TWENTY  NINE  PALMS. 
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TABLE  4 

APPROXIMATE  ANGULAR  DEFLECTIONS  (Degrees)  AS  ESTIMATED. 
From  the  Muzzle  Blast 


Event 

1 

2 

3 

4 

3 

-50 

- 

-50 

-35 

5 

-40 

-40 

-50 

-60 

7 

+40 

-40 

-40 

+15 

8 

+60 

-35 

-60 

-10 

9 

0 

-35 

-55 

+5 

10 

+60 

-40 

-50 

-10 

11 

-50 

-40 

-60 

-15 

13 

-15 

+80 

-55 

+10 

15 

-10 

+80 

-60 

+50 

16 

-10 

+80 

-60 

+60 

19 

-10 

-20 

-60 

0 

20 

-10 

-25 

-50 

-10 

21 

-10 

-30 

-50 

-10 

37 

-75 

+40 

-50 

+40 

38 

- 

+40 

-40 

0 

Range  of 

True 

Az Imuths 

-35 

-26 

-13 

1 

-52 

-42 

-36 

-26 

Average 

of  Estimates 
-8.6  1.1 

-53 

2 

5 

6 

]_ 

8 

9 

-50 

-30 

0 

-40 

-50 

-10 

-30 

-50 

+60 

+60 

+50 

0 

-30 

+60 

+40 

+50 

0 

-15 

+50 

+40 

+50 

+60 

-15 

+45 

+35 

+50 

+0 

-15 

+45 

+40 

+60 

+60 

-20 

+45 

+35 

+10 

-10 

+60 

+60 

+45 

+15 

-15 

+60 

+50 

+40 

+10 

+60 

+60 

+45 

- 

0 

+60 

+60 

+45 

+40 

-30 

+60 

-10 

+50 

+40 

-10 

+50 

+60 

+40 

+40 

-10 

-15 

-60 

+60 

+50 

+10 

+10 

+80 

+10 

+50 

15 

26 

+36 

+45 

+52 

-15 

-1 

13 

26 

35 

13 

17.3 

11 

39 

36 
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6.4  CORRELATION  ANALYSIS 

The  location  method  based  upon  relative  time  differences  between 
mortar  recoil  and  Impact  requires  a time  measurement  between  signals 
which  are  recorded  by  the  same  sensor.  Likeness  between  these  signals 
because  of  similarity  of  wave  path  suggests  use  of  cross  correlation. 

A processing  program  was  written  for  the  PDP-8  computer  which 
could  determine  the  lag  corresponding  to  maximum  correlation. 

To  provide  the  necessary  input  to  the  location  program  correlations 
were  made  between  each  explosion  event  of  the  cluster  which  surrounded 
each  mortar  firing  position.  Correlations  were  made  over  windows 
which  expand  the  entire  seismic  signature.  Each  correlation  function 
is  represented  by  128  lags  of  one  sample  each.  Since  the  largest 
signals  are  due  to  surface  waves  and  represent  the  majority  of  seismic 
energy  the  correlation  is  essentially  one  of  surface  (Rayleigh)  waves. 
The  similarities  between  signals  was  found  to  be  reflected  as 
similarity  between  correlation  functions.  Superposition  of  correlation 
functions  indicates  the  path-time  differences  to  each  sensor. 

Figure  31  shows  the  overlay  of  two  correlation  functions  for  sensor 
positions  2 through  7.  The  corresponding  differences  in  path 
length  (meters)  for  the  respective  sensors  are  -12.0,  -8.0,  -3.0, 

2.0,  7.0,  11.0,  and  15.0.  These  shots  were  separated  by  only  25 
meters.  The  step  out  corresponding  to  change  in  path  length  is 
clearly  evident  for  each  channel  and  demonstrates  inherent  capability 
of  the  system.  Figures  32  through  34  show  further  examples  of 
correlation  overlay  for  the  Twenty  Nine  Palm  data.  Figure  32  contains 
the  autocorrelation  of  shot  22  which  was  detonated  at  the  1.5 
kilometer  firing  position  and  the  cross  correlation  with  the  recoil 
seismics.  Figure  33  shows  the  relative  comparison  of  three  events 
separated  by  125  meters.  The  dash  lines  point  to  the  changes  in  align- 
ment with  corresponding  changes  in  path  length.  These  shifts  all  amount 
to  less  than  a one  cycle  change.  Figure  34  shows  a similar  case,  but 
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FIGURE  32. 


OVERLAY  OF  SEISMIC  CORRELATION  FUNCTIONS  FROM  RECOIL  AND 
EXPLOSION  EVENTS  RECORDED  AT  THE  SAIiE  POSITION.  TOP  TRACE 

» MIDDLE  EV.  38  (RECOIL) 

WITH  EV.  11.  BOTTOM  EV.  37  (RECOIL)  WITH  EV.  11.  ^ 
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FIGURE  34.  OVERLAY  OF  SEISMIC  CORRELATION  FUNCTIONS  SHOWING  CH.\NGES 
IN  ALIGNMENT  FOR  SELECTED  SENSOR  POSITIONS.  TOP  TRACE 
EV.  22  WITH  EV.  22,  MIDDLE  EV.  24  WITH  EV.  22.  BOTTOM 
EV.  23  WITH  EV.  22,  SEPARATIONS  0,  125,  50  METERS 
RESPECTIVELY. 
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with  greater  path  differences.  Note  that  the  proper  alignment  does 
not  always  Involve  the  peak  of  each  trace,  that  is  the  peak  correla- 
tion cannot  be  used  exclusively  for  purposes  of  obtaining  event  time 
difference.  Examination  of  the  computer  selected  peak  confirmed 
it  to  be  an  unreliable  method.  In  order  to  obtain  the  estimated 
relative  lag  times  between  recoil  and  impact  we  were  forced  to  make 
visual  examination  of  correlation  function  sets.  While  such  analysis 
would  not  lend  itself  easily  to  auto  processing  it  may  be  possible 
to  develop  suicable  programs. 

Seismic  correlation  functions  and  signature  data  are  given  in 
Appendix  3. 

6.5  CALCULATION  OF  LOCATION  ESTIMATES 

The  method  of  differences  algorithms  (seismlcs  on  seismics)  was 
modified  to  accept  the  processed  field  data  rather  than  simulated 
times.  The  program  continued  to  allow  for  sensor  position  and  return 
fire  errors.  Location  estimates  were  made  sequentially  but  without 
any  dependent  adaptive  structure. 

Location  estimates  were  obtained  for  each  mortar  firing  position 
using  the  individual  explosion  (calibrating)  data.  Estimates  were 
calculated  using  relative  times  based  on  acoustic  and  seismics  arrival 
picks  and  correlation  alignment  and  were  averaged  over  the  sensor 
array.  Computer  output  is  suinmarized  in  Tables  6 though  10  which 
give  event  number,  separation  distances,  and  location  error.  The 
seismic  times  from  correlation  in  general,  produced  very  poor 
results  probably  due  to  picking  the  incorrect  correlation  peak. 

Outside  of  path  velocity  variants,  the  correlation  time  differences 
should  smoothly  increase  or  decrease  from  sensor  to  sensor.  When 
analyzed,  these  rimes  were  found  to  be  often  erratic  which  could 
easily  increase  the  measurement  error.  Times  were  smoothed  to 
principal  trends  independent  of  actual  path  distances.  Location 
errors  resulting  from  rerunning  this  data  are  listed  in  Tables  11  and 
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TABLE  10 

LOCATION  ESTIMATES  FROM  SEISMIC  DATA  WITH  SMOOTHING 
TWENTY  NINE  PALMS  SITE  (1.0  km  site) 


Event 

Separation 
Distance  (m) 

Location  Estimate 
(x,y) 

Error  (m) 

1 

250 

1036,1253 

267 

3 

125 

996,963 

37.6 

5 

250 

921,1185 

201 

7 

125 

1004,961 

39.0 

8 

125 

1003,1022 

22.4 

9 

125 

1019,911 

90.7 

10 

50 

1014,1022 

26.1 

13 

125 

1031,1537 

538 

15 

50 

981,1004 

19.5 

16 

125 

706,1697 

757 

19 

125 

952,1192 

198 

20 

125 

975,1112 

114 

21 

250 

1182,632 
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TABl.E  .11 

LOCATION  ESTIMATES  FilOM  SEISMIC  DATA  WITH  SMOOTHING 
TWENTY  NINE  PALMS  SITE  (1.5  km  site) 


Event 

Separat ion 
Distance  (m) 

Location  Estimate 

(x.y)  - 

Error  (m) 

23 

64 

1075,1508 

81.7 

24 

131 

1228,1787 

33.6 

63 

253 

857,1242 

330.5 

64 

131 

952,1760 

225 

65 

103 

933,2131 

594 

66 

156 

1138,2189 

664 

67 

156 

1194,2025 

523 

68 

101 

975,1614 

78.1 

69 

258 

1277,1870 

431 

70 

165 

1168,1666 

210 

71 

290 

958,1716 

180 

73 

133 

928,1712 

188 

74 

85 

868,1399 

193 

75 

210 

979,1098 

443 

76 

327 

1234,2334 

S27 

77 

327 

842,2683 

1154 
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12.  Results  were,  improved  especiallv  at  the  1.0  km  site  at  Twcntv 
Nine  Palms. 

6.6  DlSC'JSSlOr^  OF  RLSULTS 

While  the  location  errors  presented  in  this  analysis  are  not 
encouraging,  they  help  point  out  some  of  the  problems  in  using  a 
seismic  approach.  The  basic  difficulty  is  finding  a technique  to 
make  accurate  measurements  of  the  time  difference.  Using  simulated 
sensor  position  and  return  fire  errors  did  not  greatly  affect  the 
results.  It  seems  reasoitable  tc  conclude  that  the  prime  reason  for 
large  location  errors  was  due  to  InaccuracieE  in  timing  the  seismic 
signals  and  perhaps  variations  in  seismic  velocity  within  the  near 
region  around  each  mortar  firing  position  as  defined  by  the  star  cluster 
of  explosion  sites. 

The  location  algorithm  was  simulated  as  discussed  in  section  3 under 
the  assumptions  of  constant  near  velocity  and  a timing  measurement 
error  which  was  proportional  to  the  separation  distance  between 
the  enemy  mortar  and  return  fire  impact. 

In  the  simulation  process  as  improved  location  estimates  were 
made,  the  measurement  errors  for  the  next  return  shot  were  reduced 
allowing  convergences  to  the  desired  position.  As  values  of  O'p  were 
Increased  'ihe  convergence  process  required  an  additional  number  of 
return  rounds  and  more  frequently  did  not  converge  at  all.  A value 
of  O.J,  = 10  ms.  was  found  to  be  necessary  to  have  acceptable  results 
Since  the  actual  velocities  were  about  twice  those  used  for  purposes 
of  simulation,  the  acceptable  value  of  o.^,  could  be  nearer  to  5 ms. 

Wlien  simulations  were  made  without  the  latter  assumption,  changes  of 
convergence  were,  greatly  reduced.  In  this  case  the  multiple  shot 
advantage  is  reduced  to  statistical  improvement  of  an  average  location 
estimate.  Convergence  is  by  improvement  of  an  average  rather  than  by 
error  reduction. 
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Examination  of  the  tables  of  location  errors  seems  to  indicate 
that  the  measurement  error  is  not  related  to  separation.  The 
measurement  error  which  occurs  with  real  data  is  a combination  of  both 
error  in  technique  and  error  due  to  variations  as  a result  of  the 
local  geology.  In  particular  the  signal  may  show  slight  changes  in 
frequency,  amplitude,  and  phase  velocity  which  would  contribute  to  the 
error  in  measuring  time  differences. 

Using  the  first  sei.smlc  arrival  as  a means  of  measuring  relative 
time  difference  gave  generally  poorer  results  than  using  correlation. 

The  actual  first  arrival  is  difficult  to  identify  because  of  low  signal 
levels.  The  calculated  first  arrival  travel  times  when  combined  with 
known  surface  path  lengths  did  not  yield  constant  velocities.  First 
arrival  velocities  for  the  Twenty  Nine  Palms  data  were  about  1800  m/sec 
and  exhibited  several  percent  variation  in  travel  time  to  the  same 
sensor  as  between  shots  with  corresponding  minor  changes  in  path 
length. 

The  acoustic  time  delays  produced  superior  location  estimates  to 
those  based  on  seismic  first  arrivals  mainly  because  acoustic  arrivals 
could  be  easily  identified  to  within  one  or  two  samples  (3  to  6 ms). 
Acoustic  velocities  while  affected  by  met  conditions,  are  perhaps  fairly 
uniform  at  these  short  ranges.  Acoustic  time  differences  were  six  times 
larger  than  the  seismic  and  therefore  less  sensitive  to  measurement 
error  in  determining  location  estimates. 

Cross  correlation  of  the  seismic  signature  is  considered  to  be 
better  than  using  first  arrivals  but  is  not  without  problems  and 
limitation  which  require  further  study.  The  seismic  surface  waves  were 
observed  to  have  a principal  frequency  content  between  20  and  30  hertz. 
At  greater  ranges  principal  frequencies  would  be  expected  to  be  lower 
due  to  selective  absorption  and  dispersions.  The  signals  and  their 
correlation  are  largely  sinusoidal.  For  this  reason  it  is  difficult 
to  pick  the  proper  correlation  peak  for  best  signal  alignment.  Use 
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of  the  maximum  peak  was  found  to  be  completely  unreliable  since  such 
peaks  had  often  only  slightly  greater  amplitude  than  other  lags.  A 
single  cycle  error  in  alignment  can  represent  a 30  to  50  ms  error  which 
is  unacceptable.  If  proper  alignment  could  be  achieved  a measurement 
error  of  3 to  10  ms  could  be  expected  at  a sampling  rate  of  300  samples 
per  second.  A higher  rate  would  tend  to  reduce,  this  error. 

This  improvement  in  location  accuracy  that  was  achieved  by 
smoothing  the  time  differences  for  the  1.0  kilometer  site  at  Twenty 
Nine  Palms  shows  that  some  method  is  needed  which  collectively  aligns 
all  sensor  signals  simulatneously . The  poor  results  achieved 
otherwise  indicates  that  correlation  by  itself  is  not  adequate.  While 
it  is  now  known  why  the  same  smoothing  process  did  not  improve  results 
at  the  1.5  kilometer  site  it  is  recognized  that  the  time  differences 
were  much  less  and  therefore  were  sensitive  to  alignment  errors.  In 
this  case,  a greater  sensor  spacing  (base)  may  be  required  to  achieve 
the  same  degree  of  location  accuracy. 

It  should  also  be  mentioned  that  location  estimates  were  based 
upon  an  average  of  all  possible  sensor  solutions  and  therefore  if 
certain  time  differences  were  in  great  error  they  would  tend  to  throw 
off  the  entire  estimate.  Some  improvement  might  result  if  a weighted 
geometric  average  rather  than  an  arithematic  average  were  used  to  make  the 
estimate.  Missing  sensor  data  are  indicated  in  the  tables  of  Appendix  B. 

While  the  field  data  was  not  processed  to  be  used  with  the 
hyperbolic  location  algorithm,  the  resux*'s  of  using  this  location 
techniques  would  not  likely  be  an  improvemerii;  over  those  estimates 
reported  above.  The  basis  for  this  conclusion  is  mainly  due  to  the 
fact  that  greater  correlation  of  seismics  was  observed  between  shots 
than  between  sensors. 

In  summary,  it  is  felt  that  significant  advances  must  be  made  in 
the  methods  of  observing  seismic  time  delays  before  the  advantages  of  a 
seismic  location  system  can  be  realized. 
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7 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  and  recommendations  are  made  after 
completing  this  study  of  seismic  mortar  location. 

7.1  GENERAL  CONCLUSION 

While  this  study  found  some  promising  characteristics  of  a 
seismic  mortar  location  concept  sufficient  limitations  at  this 
point  of  development  will  require  further  exploratory  development 
to  resolve  before  beginning  any  advance  development  work. 

7.2  SPECIFIC  CONCLUSIONS 

Range  Limitations  - Under  quiet  background  conditions  at 
Twenty  Nine  Palms  seismic  signals  with  low  noise  levels  were  observed 
from  81  miri  mortar  recoil  to  a range  of  about  1800  meters.  Events 
at  greater  ranges  were  planned  for  the  field  measurements  program 
but  were  not  completed  because  of  time  and  windy  conditions.  It 
is  felt  that  seismic  signals  from  recoil  could  be  observed  to  ranges 
of  3000  meters  under  quiet  conditions.  The  use  of  signal  enhancement 
techniques  could  double  this  range.  At  Camp  Pendleton  under  high 
background  noise  conditions  generated  by  base  activities  reliable 
seismic  signals  from  mortar  recoil  could  not  be  observed  at  500  meters 
range.  While  loss  of  signal  could  be  due  in  part  to  path  effects 
and  energy  coupling  problems,  low  signal  levels  were  likely  due  to 
the  high  noise  levels.  Required  field  gains  were  20  dB  greater 
at  Camp  Pendleton  than  at  Twenty  Nine  Palms. 

System  Design  Problems  - Although  only  a cursorary  examination 
of  system  design  was  made  no  serious  problems  could  be  identified  which 
would  not  allow  building  a system  with  off  the  shelf  technology. 
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Basic  Information  Requirements  - Tha  essential  :Input  information 
needed  to  establish  mortar  locations  using  selsmics  Includes 
positions,  impact  positions  of  directed  fire,  and  the  seismic 
signature.  Of  primary  importance  is  the  use  of  directed  fire  for 
calibration  and  ranging.  Thus,  this  study  could  not  find  any  practical 
means  of  obtaining  necessary  calibration  information  covertly. 

Signal  Characteristics  - Seismic  surface  waves  (Rayleigh  waves) 
as  observed  on  a vertical  axis  seismometer  were  found  to  have  the 
largest  amplitudes.  These  surface  wave  signals  were  identical  for 
recoil  and  explosion  events.  Also  shots  fired  in  the  vicinity  of  a 
mortar  position  showed  a great  deal  of  similarity  as  did  the 
corresponding  cross  correlation  functions. 

Algorithms  - Of  the  two  algorithms  investigated  the  method  of 
relative  differences  (seismics  on  seismics)  appears  to  be  the  most 
promising  for  a small  light-weight  system  concept. 

Location  Accuracies  - Location  accuracies  produced  by  utiliza- 
tion of  field  data  were  generally  disappointing  for  the  seismic  data. 
The  acoustic  data  appears  at  the  short  ranges  under  consideration  to 
produce  superior  location  estimates.  Such  differences  could  be 
attributable  to  the  fact  that  the  seismic  velocities  were  five  to  six 
times  faster  than  the  velocity  of  sound.  The  higher  velocity 
media  implies  a greater  sensitivity  to  timing  errors.  Some 
difficulty  was  found  in  timing  the  seismic  waves.  The  use  of  peak 
correlation  was  found  to  be  an  unreliable  means.  Location  errors 
obtained  irom  timing  seismic  waves  were  not  related  to  the  separation 
distance  Leuv,^cen  a known  impact  and  location  of  an  enemy  recoil.  The 
timing  error  was  apparently  not  related  directly  to  the  separation 
distance  between  recoils  and  impacts  as  expected. 
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The  location  estimates  made  using  acoustical  and  seismic  data 
recorded  at  Camp  Pendleton  were  not  useful.  The  poor  quality  of 
these  estimates  may  be  attributable  in  part  to  the  high  noise 
levels  and  frequent  interferring  sources  from  base  activities. 

7.3  GENERAL  RECOMMENDATION 

It  is  recommended  that  further  exploratory  work  be  considered 
for  development  of  a light  weight  man  packed  seismic  system. 

The  system  concept  should  be  expanded  to  targets  other  than  mortars 
and  should  fully  exploit  the  use  of  acoustic  signals  as  well  as  Uie 
seismic . 

7.4  SPECIFIC  RECOMMENDATIONS 

1.  Investigate  and  evaluate  using  existing  data  methods  of 
signal  enhancement  which  would  be  applicable  to  this  system 
concept.  Project  maximum  effective  range  for  various 
targets . 

2.  Determine  surveillance  capabilities  of  the  seismic-acoustic 
location  system  against  such  battlefield  targets  as  trucks, 
tanks  and  a walking  man. 

3.  Investigate  and  evaluate  means  of  signal  processing  which 
will  produce  accurate  signal  timing  information. 

4.  Investigate  methods  of  improving  the  signal  to  noise  of  the 
saismics  by  analyzing  the  collected  data  in  the  frequency 
domain . 

5.  Investigate  the  use  of  acoustics  for  location  as  well  as 
the  advantages  of  seismic  for  identification  or  discrimina- 
tion of  a mortar  recoil  from  a shell  impact. 
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APPENDIX  A 

FORTRAN  Listing  of  Simulation  Algorithms: 

1.  Hyperbolic  Location  with  Kalman  filter. 

2.  Subprograms  to  the  Hyperbolic  Location  Algorithm. 

3.  Location  Algorithm  using  the  Selsmics  or  Seismic  technique 
(subroutine  k4). 

4.  Histogram;  Same  location  simulation  as  (3)  but  designed 
to  make  many  iterations  and  printout  a histogram  of 
location  errors. 


97 


,n  n o 


FORMERLY  WILLOW  RUN  LABORATORIES.  THE  UNIVERSITY  OF  MICHIGAN 


c 

C PROGP**!  TO  LOCATE  TJIPGFT  HORTAR  FIRE  POSITION  IMFOSPH  A SncCFSTlOH  OF 
r APPE0TT»1AT10NS  HTILIZINS  TARGET  PTRF  AND  N SF.NSOR  LOCATIONS 


C 

C V 

C 

C TPW 

r X (T)  ,Y(T) 
C XA,YA 
C yp,YP 
C XR,YR 
C SG 
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C SP 


ARIAtJLFS  ttSFD  BY  FFOGRAM 

SWITCH-:  foh  fixed  SFNS0FS,  = I FOn  rlRF  SENSORS  AUOrJND  POINT 

ACTUAL  S'’NSOB  LOCA'^IOK  COOPDT>!  ^l[>F,S  PFLATIVE  TO  CO'IPUrEP  POSI"'10N 
APPPOXIFAT’ION  OF  FTKXnR  P0IN"T 
ACTTJAL  FIRING  LOCATION 
BETBRN  FIRE  SOURCE 

STA'IDAPD  DEVIATION  FOR  FPRO'>S  IN  DI3T.AVCF  M ^ASU H FH FM  T I 

STD  DFV  FOR  TINE  F’^ASUR EN^NTS 

STD  EFT  FOP  DEPIFCTION  IN  PF'’UP;i  FIR” 
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<^ALL  Tr"E  (7,  “,INrT)  ' 
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CAIL  GPAVDI  (iriTT) 

r TN'r!!"  VEL0CTT7PS  AND  S'^AUDAF^  VIA '"ION' 

C 
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TF  (V.GT.M)  GOTO  9^'' 
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Y (4)=Y(2) 
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Y (9)  =y.*--SPTL 

X (7)  =XA 
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Y(P)=Y{6)  ___ 
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7*='  (PHT.L".-  .)  pni=ri+r»<i 

pp=p»-  /2 . -f”'T 

yr  (TJ  ="'  ftri  V (P'!I)  (PF)  +X  ( I) 

u yp(T)  f7n*ccs(o^)  ♦•y  (T) 

'',r>"’0  6 


•’"A?  ((X (I)  ,y 

XF  F’^An  (0:{T)  ,Y  (I) ) ,T  = *. 

DO  I**  1=1,^' 

vr«  (X)  =r.r  M;n  (P*:,  X (T)  ) 

1'  yp  (T)  =o<?n'"' (SP,  Y (I) ) 

PFT  iSP  C0N7TA*”‘P 


N> 


•1  = v-t 
ro'-  = 1 . z"’  - 
or  in  1-1,'’' 

1.’  o(7,T)-''r‘i 


r*?\'Fo|iT-.  a’‘-*’AY  ('F  TAPO-T  rnp  TMT'^TAI,.  P7"1JFM  FIPi: 


o ''  .1  r 7 ( 1 ) =Y  A 

VT  (1)  =vs 

no  11'“  T=1,n 

XTC'  + I)  ~V'* 

XT  ([♦  1)  --y\ 

XT(t+7)  r XT  (•?) 

V"  (T*.  2)  = Y."  ♦■3TF? 

XT(T*3)  =XA 
YT  {'*^)  =Y--^  (3) 

Y7  (T*U)  =Xt-PTT’’ 

VT  (T*  u)  - Y"  (■*) 

XT  (T*  R)  =7T  (F) 

Y""  (I  + 'i)  -YA 
YTfT+r.)  rrXT(T) 

YT  (1  + 6)  -YA-PTFP 
XT  (x*7)  -YA 

YT  (T  + 71=YT(7) 

XT  (X  ♦«)  ■=  X"'  (?) 

YT  (T+H)  = YT  (X) 
sxt'p*  n’-'L 

IF  (T?ir.  LT.  I'")  ''•OTO  ^2" 

IK  rONTT  «]u: 
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120  COVTTUnF 

c 

C CALr  nPLTfi  ’^'S  F0!>  iNTTTftl  KFCOTL 
C 

W(l)  =D(X(1)  ,Y(1)  ,Kf  ,Yn 
T (1)=W{1)/V?  (1) 
no  2''  I = 2,'J 

W(T)=D(X(T)  ,Y(T)  ,XF,Vn 
T(T)  =H(T)/VP  (T) 

20  0T(T)='^  (T-1) -T(T) 

m(1)  =T(N)-”(1) 

DO  13  T=1,M 

V(l)  =1./V1  (T) 

13  VP  (I)  =V1  (I) 

K=" 

WPIT?  XF,YF,X!l,YA,SS,ST,SR,S:^ 

BPTTP  (d,X(T),Y(T),XE(M  ,YE{I)  ,VR(I)  ,V1  (T)  ,DT(T)),T  = 1,N) 

no  21  1=1, M 

21  DTE(T)=D*(T1  ♦OFAMD(nT,<‘.) 

WFITP  (6,25")  ( (I,DTE  (D)  ,I  = 1,N) 

r 

C IMPROVE  VELOCT"'Y  ES'^IMATFS  "‘RmtlRH  MULTIPI.T  T'»PAC'^S 
C 

T.=^  X 

DO  111  1=2, IMP, 2 

PR=0nA?!D(SP,''.) 

PD=G3>HD  (S0,0.) 

PHT  = ATAN  ( (XT  (I)  -XR)  / (YT  (T)  -YR) ) 

IP  (PHT.LT.C.)  «>HI  = PHI*PI 
PP=Pl/2. -PII 

Xl=Pp*STN  (PHI)  ♦PP*SIN  (PP)  »XT  (I) 

(Re:>)«‘Yi?.  (IJ, 

C - 

Ifi'*  PR=GRANn  (SP,':  .) 

PD=r.R  AND  (S0,n.) 

PHT=ATAN  ((XT  (1-1)  -X R)  / (YT  (I- 1)  -YR)  ) 

IF  (PHT.LT.O.)  nHi=PHT+PI 
PP=PI/2.-PHI 

XTI  = PR*SIN  (PHI)  ♦Fn*STN(PP)+XT(T-1) 

YIT=PR*CO.S  (PHI)  +FD*C03  (PP)  +YT  (1-1) 

IP  (XI.  FO.XII.AKD.YT.EO.YIT)  f',0’>'0  ISO 

FIND  TR.AVET.  TI'^pS  AND  VELOCITIES  FEOr  EACH  I.MPACT  TO  F.ACH  SENSOR 

DI(1)=D  (X(1)  ,Y(1)  ,XI,YI) 

0P(1)  =n(X(i)  ,Y(1)  ,xir,Yii) 


Reproduced  From 
Best  Available  Copy 
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no  1U'’  .1^^/ 

OT  (J)  =r.  (X  (.7)  ,v  ( T) 


:i  ,y^) 


,YTT) 


J^P(.7)  (X  (•>)  (.7) 

■7  f ;j_  1)  =ni  (.7-  7)  ,/Vl-  (.7-1  ) -PT  (.7)  /v>  (.7)  » !•. 

7 =rn  /VQ  (.7-1)  -n?  (.7)  /V'  (I)  + 

n (j-1  ,.7-i)--^c^m.7-i)  ,Y-^  (.7-1)  ^x  Mn  .Y-n)) 


AfP("7','  .) 
:.1Nn  (ST,‘  .) 


" {.7-1,.M  = 

H ('*♦.7  -1  , .7 
" (”♦.7-1  ,J)  =-P  (X 
lU'’  ro‘iri>'iiP 


(xr(.l)  ,Y  - (.7)  ,>T(t)  ,v-(-)  ) 
1)  =•>('' <^  (.7-1)  ,Yr  (.7-7) 


) ,Y : (T- 1)) 


(.7)  rXT(T-l)  ,'^’(7-1)  ) 


CALT.  K;L1.\V  C!,V  , , P , 9 , ^1,  M) 

7.  =1.  * ?. 

•7C  13''  .i=i,r 

13'  VF(.7)  =1. /V  (.7) 

ypT--?  7. 

HPTT-  (1^,1'^')  ( (.7,V-*  (J)  ) ,.7^1  ,N> 

131  c 01“^  1*171=: 


r 

C 

r 


C 

r 

r 

C 


r 


ij*;"  gcT'TS'v  V ■;7.0''f" *■  *'1  A ‘1 5 P-^CHL  T'l=‘7>  T'lP  ‘1 7’.'  POTiJi’  7Srir’.\'’”^ 

3*i  K=K*  1 >. 

Tr  (K  .OT.  K***.  V)  '7-7TO 


r!?"  A prpoyT  vaTTP"!  ti”  '=\!7n?T  rO'>T'=T7.7‘J 

USTHG  EACH  PAT?  <1’’  OOTMT-S  Tl*:!  SSTIIA'^F  FOH  AV“?Ar,I*JS 

T.='  ^ 

yn  = '“ . 

Yn=  ' . 

nn  ■»''  Is.?," 

C.M7.  HLOCCPT-  (T-1)  ,Y7='(I-1)  ,V’-  (1-7)  , XA  , YA  , XL  , Y L , C 3 ' ) 

Xn  = v'1  ♦VT. 

Yr-Y"*YL 

L=l.*1 

3"  CCrTTN'I? 


2 SPECIAL  CAS^"' 


AA (1)  =y-  (V-7) 
\A  (?)  (*0 

AA(3)  =:<^  (1) 
Tn(1)  =Y9  (V-  7) 
m (?)  =Y”  (*•) 

I\T  (1)  =VT^  (1) 

rc(1)  =V”  (M-1) 

CC  (2)  =V=-  (V) 


«fe6/ec'“'" 
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CC  (^)  -VF  (1) 
rD(1)  =0-^^  (T) 

OP  (2)  =0T''  (1) 

r ALL  HLCC  ( A A , 00  ,rc  , 00  , XA  , T'  , I. , VT. , >*.  3 1 ) 

xn=xo+xi. 

Yn=YP  4-YL 
L=I.+  1 
r 

3^.  a\(l)=yr(s} 

t A (2»  =XF  (1) 

AA(T)  (?) 

pi?(1)=yr(v) 

= 0 (?)  =y^  (1) 

f«3(3)  =YF  (2) 
rC(1)  =V-  (^f) 

CC(2)  (1) 

rc  (3)  =Vr  (2) 

PD(1)  =DTF  (1) 

no(?)  =0'"’^  (2) 

CALL  HLOC  (A.  A,nP,CC,PD,XA,y>  ,XL,yL,63.1) 

X3=T3*YT.  ^ 

YP=Y3+yi. 

L=T  *■  1 

32  IF  (L.VO.')  -30 TO  93' 

X.\=TB/L 
Y\=Y« A 

r 

r OTFor.Ay  «em  T.-nTTf<:A""K 
r 

niJF  = P (XA  ,YA  AF,Y^) 

(f^,?''-')  K,T  ,XS  ,YA,P03 
C 

c ^A'iro?!iy  rrcTO”  r oata  ok  get  r'^u’Acr  oa^a 

c 

nfiy.=npAVDii. ) 

TF  (Onr.GT.  ■'.)  "JOTO  36 

r 

C fJE'IFTA'^F  T'^^ACT  DATA  FOP  (XA,Y*-)-GEt  poj'JF 

C 

ftC  Y=K  + 1 

IF  (Y.O^-.KTAX)  ^O-O  9:' 

F'’=GPAN'0  (F'’,T  .) 

rn=ar  Avn  (Sn,  '>.) 

^ f'T  = A TA  N ( ( X?  - XP ) / (Y  A - Y*» ) ) 

(DHT.  L'T’.'*  .)  ?FI=PHT4-PT 

''P=PT/2.  -p"T 


'S'SSs;'. 


^^flabu 


* 
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XT=  rR«5;rN  (P!'  f)  ♦•"P*?:  IM  (PP)  *X  \ 

YT^F’^^ros  (■■''•T)  +pp*co5  (pp) 

C. 

c T'-'P^rT  POP  (X.A  t PF  I,,  V Af  PP  I.) 

C 

Tflr- V A +PF  T. 

YP,.-yA  s-r^i;  I 

f,*’  rr  = GPAnn  ^ 

TPr-GO  A*'P  (F^,  “ . ) 

PHI=,^"•A^’  ((<P-XD/(Yn-Y?>)) 
tP  (PPI.T.*".  .)  Piir=PMl4-PT 

r>n=o  ▼ _n 

YTXsr-Ti  + s rU  (P!'I)  ♦■PP+F  TN  (PP)  +X'^ 
y TT  = PP*CO«;  (P»i  T)  ?T  (FP)  +Y  P 

Tf  fxj.rn.YTT.Avn.Yi.Fo.YiT)  cnri  ' 

FTVP  ""rAYFI.  Ai^ri  V F.r.OC Fl-'f'’'.  “AC:f  Tf-IPACT  TO  SENf?.')!’ 

PT(1)  =r  (XC)  ,Y{1)  ,XT,Y-) 

np(i)  =n(-f  (1)  ,y  (i)  ,xtt,yti) 
no  £ip  n 

PI(T)  (X(n  ,Y(T)  ,XI,YT) 

PP(n  =0{X(')  ,Y  (T)  ,X  JT,YTT) 

T (T-1)  =nt  (r- 1) /vr  (t-i)-pt  d)  ♦';'a.'iT  (sr, ' .) 

Y (F  + T-l)  =DP  (T-1)/VR  (1-1)  - PP  (T)  /VR  (T)  +onA‘!D  (ST,'  . ) 

(T-1  ,T-  1)  =P  (•<'■  (T-1)  ,yr  (T-1)  ,YA  ,y;,) 

(T-1  ,T)  =.n(YT(7)  ,y-r  (I)  ,y.  ,ya) 

•F  + = r (yr  (X_-I)  ^YT  (T-  1)  ,xn  ,Y'^) 

F'  (l*'*T-1,T)  --P  (.yp  (T)  ,Y^  (T) 
un  CON'*'T  PFF" 

TAYL  KAI.  (»’ , V , F. , P , ?(  , 2**1 , N') 
no  S'"  T=i,v 

t;.''  VP(T)  =1./V  (T) 

WPIT”  (6,V') 

vnjTT^  ( (i,vy:(I)  ) ,I=1,N) 

r 

r GKT  NYW  POINT  p.'^TIYA':^^ 

n 

T — * 

xn=r-. 

YM='' . 

PO  70  T = .7,  -* 

C.AT.L  FFTOCCvr  (T-1)  ^YTd-l)  ,yP(I-1)  ,*lTF(r)  , XA  , Y A , XL  , Y I, , P.?-) 

YR=xn  *yT  » • / 

ve=YR+YL 

I.=L*1 

7'“  ccnttnij" 


iOA  Reproduced  From 
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C 2 SP^CTAl,  CASFS; 

r 

A A (1)  *XF  (B-1) 

A A (2)  *XE  (N) 

AA  (.A)  =XF(1) 

BB(1)  =YP  (V.~^) 

BB  (2)  =YF  (N) 

BP  (3)  =YF  (1) 
rC(H  =VF  (N-1) 

CC  (2)  *VF  (N) 

CC  (3)  =VP  (1) 
nc(  1)  =DTF  (N) 

nn(2)  =0'^’!’  (1) 

CALI  BLOC (AA,BB,CC,nn ,XA, YA ,YL,YL,S7 1) 

XP=Xn*XL 

Yr=Yn+Yi. 

T,=I4-1 

V 

71  j\  A (1)  (N) 

AA(2)  ^XF  (1) 

AA(’)=X’^(2) 

BB  (1)  =YF  (N) 

BB  (21  =Y"  (1) 

BB(?)  =YF  (21 
CC  (1)  =VF  (M) 

CC  (2)  =VF  (1) 

CC  (31  (21 

nn  (11  =0-"^  (1) 

nD(2)  =nTw  (7) 

:AI  L lU.OC  (AA,PB,CC,Pr.  ,XA,  YA,XL,  YL,S7  ?) 
XB-XB+XL 
VB=YB  +YL 
L = T,+  1 

72  IF  (1,  .F0.''1  ^OTO  93 1 
X A=XB/T, 

VA=YB/L 

C DI?.?'^T,AY  NFW  APPPOX IMATION 
m)X--r  (XA  ,YA  , XF,  YF) 

WPITF  (P,?1F)  L,XA,YA,DnM 

CO'^0  25 

9 O'*  HPTTF  ((:,5"'^) 

ocq  CONTI  in? 

FTOP  P 
9P"  STOP  1 
9P1  STOP  2 


')  ^OTO  93  1 
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rSr!  oei-ta  t (hlc.il);/) 

PO-'^KOP^A?*'!*.*'  •.n,'  rproil.  DATA  — NSW  R!?T1NATK  BAf^ED  ON  *,T3 

SFHSCN  5?PTS=  ' , ’^’ci.  1 , * r ' »FB.  1 » 

t iriS'IR''CP  "TO  F=  * » FP  . 2) 

P1“  VOFFAT  ('C  PPW  PSTINATT:  nA5?n  ON  »,Tl 
■?,'  SENPOR  SE1’S=  ' , . 1 » ' » * » 
f, 'PIPIAHCE  TO  r.OtlRCE= ',F8  .2) 

2«^/'  POFNAT  RPrOIL  *'FI.TA  H/  pRROR:* 

• TFPRCT  PATA  — NSW  PATt!  VELOCI  riER= ') 

3'5''  FOFFAT  ( ?'=X , H , F 1 2.  U) 

jjn’'  POFNRT  (IOP'1,2) 

U'i'i  ’'OFKAT  (3T3) 

EOF  FOPWAT  ( ' O ETOP  0*) 

To?  pS;?rt  r’  r»;TT‘.i'r”s^or.p  ..n,.  r.P-.cp.  ,n.r.«  »E»  nr.- 

J,«  VFl.orTTTES  = * ) 

END 
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f"* 

r 


0' 

.-'O  T.OCA^P  ?0T\'^  n:;”'!!;  "y'’’^'’POT.TC  FPS’CT'f'OM.'-':  riir’l 
'iP"'HnD  S’^’'?PFPT  DF'^CF>-‘^ 


PMr  POrri  MLOC  (X'  ,Y  , V , , X A ,?  A , X L,  V T. , * ) 

•5r,'  j^<ru  V 

r.Tr'-rNt.r;rnM  X O)  ,'{  ('3)  ,V  (3)  ,-{2)  ,?P  (3)  ^’^•X  (2)  (2)  ,P  {2) 

7^-YA 
Yo  = YA 

Vr  ^ 

* — I * 

•='A''''"K  = 2. 

3'  DP(1)  =D(X{1)  ,v  (1)  ,XP,YP) 
np(?)  =r.  (X(2)  ,Y  (2)  ,XP,  YP) 

OP  (3)  =P  (X  (3)  (1)  ,XP,  vp) 

no  1"  ■'=1,2 

'1  = V(T)/(V  C"+1)*PP  (T+D) 

”x  (I ) = (X'’- Y (T) ) /n?(T)  (x?-x  (1  ♦ 1) ) 

’^Y  (n  = (Y  P-v  (I)  ) /nr  (T)  -i:;*  (YP-Y  (■'■  + 1 ) ) 

1-  (J)  =np  ( r) -V  (T)  *r,p  ( r^-i) /V  (T+ 1) -T<I)  *V  (I) 

3X=?.  *F  ( 1)  (13  f 2.  i'P  (2)  *FX  (2) 

'AV=2.  (1)  ’^‘FY  (1)  +2.  (2)*FY  (2) 

'^F=D  (nx,  ly,','. 

•(L=XP-F»3X/nF 
YL=YP-K*r,Y/r'F 
;p  = F ( 1)  «F  (1)  +F  (2)  *P  (2) 

'i"  ’■  (1)  = rc'  (1)  ,Y  (1)  ,XL,YT0  -V  (l)xcp  (V  (2)  ,y  (■?,)  ,y.I,,YL)  /V(2)  -T(1)  *=/  (1) 

(2)=n  (X  (2.)  ,Y  (2)  ,XI,,YL)  -V  (.?)  (X  (3)  ,Y  cn  ,XL,YL)  /vn)  - ■!'(?)  *V  (2) 

OL=F  { 1)  ’i‘P  (1)  +P  {>)  (?) 

T.  F (fIL.  'IP)  nOTO 

PALI.  FTFP  (XP,YP,XL,  YT,,XN,Y'^,F) 

XP=XL 
YP=YL 
XL=XM 
YI,  = YN 

■TPv-I 


y = r?\  r "’!f*  K 
LOTO  UP 
n <<  =F  /2  . 

FAP'PK  = 1 . 

TF  (K.G'^..0l)  r-n':’o  3^ 
TF  (GL.GT.3.)  3”' 

P'  npTnp^i  1 

"■NO 


Reproduced  From 
Best  Available  Copy 
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IWD  E*«n  FY^INT  OF  OIKFrTPn  VFCT'M*  F *>n*«  fX1,Y1)  THROUGH  (X2,Y2) 
WHICH  IS  niSTXMCF  K FROS  (X2,Y2) 

SOBROMTIRF  *?T«fO  (Xl  , Yl  ,X2  , Y?  ,X  3 , Y1 , K » 

THPT.ICIT  (A-7.) 

TF  (X1.’'0.Y2)  00*^0 
«»(Y2-Y1)  ./(X2-X1) 

n«ROnT (K*K/ » 

X3»X2*0 

IF  (n*(>2-X1)  .IT.  ?.  > X3*X2-n‘ 

V 3=p*n*Y? 

RETURN 
10  X3=X2 
XK=K 

Y3=T2+STGN  (XK,  (Y2-Y1)  > 

RE'^unn 

RNO 


USE  KM.»'M4  FTLTFP  "O  PTND  FEW  w?‘T.OCITTrS 


SUBPOIITTHF  KM.MAN  (ll,V,Z,P,R,  1,  N) 

CC**EON  PTH’/ (2'>  , 2'') 

PFAL*fi  01, 02 

RFAL*«»  K 

DTFFNFnH  M (2'',2C)  ,T’»(2P)  ,V  (2'')  ,%  (22)  ,»>(20 ,2l  ) , R (2!  , 20)  , K (20,  2C) 
S,HT(20,?^)  ,F1  (20,2~)  ,W2  (2^,2'')  ,Dl  (2*  ,2")  ,1)2  (20,20) 

CAl.T.  ’’■RAHS  (!' ,MT  , h,N) 

CALL  f'"LT(H’’,",R1  ,M,S,F) 

CALL  MUI,'’’ {'s'1  , H,H?,H,1,N) 
no  1?  T = 1,*l 

0 n 1 .T=  1 , V 

HTU  V (T,J)  =”"NV  (’^,.1)  ♦H2  (T,,1) 

01  (T,.l)  =PIVV  (I,vT) 

CALL  urn  T(H,  V,».'2,H,N,  1) 

CALL  I*IV  (N,2'' ,o1,TP,20,r2) 
no  2"  I=1,M 
no  2"'  .1=1,?; 

20  p (T,.7)  =07  (^,.1) 

00  .30  T=1,'1 

■70  W2  (I,  1)  =7.  ("!) -”2  (T  ,1) 

CAM.  1MT.t(p,W1,K,V,*J,!1) 

CALL  rULT  (K,  W2,?H  ,N,f!,  1) 

PC  a''  T=1,*l 
«o  V (T)  = V (I)  +1J1  (T,  1) 

npTnp  n 


108 


. II 


2»i 


FOAMERUV  WILLOW  RUN  LABORATORIES.  THE  UNIVERSITY  OF  MICHIGAN 


ciinRonTTfjf:  multipi  y .-ia 

SUFPO vni,'^  (A,n»c,L,*',N) 

A (2'',2'')  ,n(2',r‘)  ,C(2' , 2 ) 
on  1'"  T-1,L- 
DC  1 ''  J=  1 , >’ 

no  " K-  1, 

??=A  (T,K)  *n  (K,.T) 

1'  C (T  , J)  = C (T,  J)  +"? 

!?FTgt?  y 


r SHBROriTT  fjy  TO  '^  r>A»IBPOS~  *1ATRrcn3 

r 

<^riR=OIJTT  n?  7FAKB  (A,P,m,N) 
nTr^wFIOM  A {2'‘'' r B ( 2'' , 2 
DO  1'"  1=1,  K ^ 

no  in  J=1,M 

FFTOPN 

FND 

r* 

C S tIFfon?T»!F  '"0  oF'^UP  ARnAYS  FOP  KAL'IAM  ''ILTFF  OF  RFCOIL  DATA 
C 

FanponTTN-F  kpse"^  (x,  y,m,xa  ,ya  ,xn,Y»,  n') 

DTFFNSIOM  X (10)  ,Y  (1‘>  ,K(1  ^ 

M =y  - 1 

DO  2":  T=1,  Y 

A1  = E (X(T)  ,Y  (I)  ,XA,YA) 

A2=D  (X  (1  + 1)  rY  (T+1)  ,XA  ,YA) 

X11=  (XA-X  (T)  ) /A1 
X2r  = (XA-  X (T  + 1)  ) /A  2 
V11=  (YA-Y  (I)  )/Al 
Y22=  (YA-  Y(T+1))  /A  2 
M (J  , T)  = y 11*  (Xn-XA)  + Y1  1Y  (YP-YA)  + A1 
2'’  w (I, T + 1)  =X22*  (XA-Xn>  +Y22*  (YA-YD)  -A2 
•■'E'rriK  N 
PND 


r.  IIAMPY  Dir^TAYCF  FUNCTION 

r 

^IINCTTON  D(A,D,W,Z) 

n=F0P7.  ( (A-F)  * (A-K)  t [2-7.)  * (P-7.)  ) 

RFFtlP  N 

Fwn 
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siibrouttnii:  k« 

0(A»B»W,2)a^QRT(  (*.W)w(  A-W)4.(R-7)*  (h-7)  ) 

HFAL  K 

OTMFNSinM  xnO),Yfin),Xlf(tO),Venn),TR{tO),TIflO),r.TnO^ 
*rS(10i»C(10i,SM(10),X?(?),V?(?),Kfa),Vlf2),YJf21rn(10»ai 

DATA  PI/3.iai'5926‘53S/ 

READ  (5, IDO)  M,IKMAX 

RFAD  rSaSO)  (Of(T),Yn)),lsl,5) 

M*N«1 

IF  (N,GT.S)  RFAD  (5, 130)  C ( Y (T ) , Y f I ) ) / 1*6, N ) 

RFAO  (5,130)  3S,ST, AR,3D, YA, Y4,yF,YF,VR, VF,Dei  ,B0 
IKaO 

INITIAlTZF  MORHAUZEO  random  NUMBER  GFNFRATOR 

call  T1me(7,0,IN1T) 

IK)JT«IMTT*2+1 

call  RRANDKIMIT) 

GFT  RFcniL  times 

nn  to  1*1, N 

T«(n»0(XfI),Yn),XF,YF)/VP 

XP(T)aX(I)TGRAND(SS,0,) 

YF(n«Yf  I )tr,RANn(RS,0.) 

10  CONTINUE 
MLaO 

RFTI.IRW  PIPE  TO  fXA,YA) 

50  IXaTK+1 

IF  (IK.GT.IKMAX)  GOTO  90 
ISsO 

XTsXA+OFL 
YTsYA+DFL 
ERbGRANO(SR,0.  ) 
eDaGRAND(SD,0.) 

PHIaATAM(XT/YT) 

IF  (PHI.'lT.O.)  PHIsPHT  + PI 
PPaPI/2.-PHr 

XIaER*8lN(PHI)+EO*8TN(PP)+XT 

YTaFR*C03(PHI)teOACOS(P»)+YT 

OUM1rO(XF,VF,XI,YT) 

XRaO. 

YRsO. 
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56 
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70 

71 
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73 
7a 

75 

76 

77 

78 

79 

60 
51 
82 
63 
8a 
65 
65 
87 
68 
69 
90 


C GPT  IMPACT  TRAVFL  TIMPS 
C 

on  30  J»»1,N 

TT{J)sD'XfJ)»TlJ)*XT,Yn/VR 
SM(J)a(YT-y6(J)  WrxT-VEf  J)  ) 

THaATANfSM(.J)) 

IP  (TH.l.T.O,)  THsTH  + PT 
Cf J)sr06(TH7 
Sf J)aSlM(THi 
SM(J)a-l,/SH(J) 

30  cnNTlNUP 

KaisTifn-TRU) 

KfaiaTI C2)-TRf2i 
DO  20  I=3»N 
on  20  Lsl»2 

TsTICn-TR{T)+DUMlB.r5RANOC9TBO.)/RO  + Ka) 

20  an  ,L)aYT-T*VF*S(T)-SM(T)*XT  + SMf  n*T*VE*CCn 
on  50  Iaa»H 
OP  70  L*l»2 

X1(L)3(R(T,L)-Bf  I-l  ,Li)  / (3M(l-n-RMtn  ) 

Yi  (U=SM(T-n*Xl  CL)+Bf  I-l  ,L) 
X2(L)aC6(Tfl,L)«B(I,L))/(8M(n-RMfI  + l)) 
Y?a)sSM(T)*X2(UfB(I,L) 

70  COnTINUF 

SMla(YH2T“Yl(n)/(yi(21-Xl(n) 
SMIIa(Y2(2)-Y2(n)/(X2(2)-x2(ni 
XTEMarvl ( t )-SHI*xl { 1) • Y2 ( 1 ) +SM J T * X2 U ) ) / ( SM T I -SM 1> 
YTEMa8MT*XTFM*Yl  (1  )-SMI*Xl(n 
xaryBAX  TE5 
50  YPsYB+'iTEM 
X6aXB/(M-3) 

YRaYB/(M-3) 

XAa(XA*(IK-t )+XR)/IK 
YAr(YA*nx-l)+YB)/IK 
DUM2rO{XF,YF,XA, YA) 

WOITE  (5,R0O)  I«,XA, YA,nUMl ,0MM? 

GnjO  SO 
90  RFTIIRM 
100  format  (2T3) 

ISO  format  ri2Fio.a) 

200  format  (I5»*  ESTs' ,F6.1 b ' , ' »F8,l. ' bDIST  TO  IMPACTa',F8 
7, ' , OTST  to  EST*',F6.1) 

END 
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SUr^POUTTv  E KL5T 

n {A  ,B  /<,Z)  =?onT  ( (3-  M)  * {A-W)  f (H-Z)  *(B-Z)  ) 

PFAI.  V 

DT''F^  0T2M  Y (1T  ) , ? (1  ''1  ,XF  (1  ")  , Y 5 (1  ( 1 '' ) 1 T ) 

r,  5(1*  ) ,r:  (r ) (1  >)  ,y  2 (2)  ,Y?  (2)  ,K  (2)  , A i (2)  , y 1 T ,2) 

F , THTfjT  (2  6) 

DATA  P-!  / 3.  1 41  5'^26535/ 

R"An  (5,K') 

’”!^An  (5,15*)  ((X(I)  ,Y  (T))  ,T-1,5) 

1 = M-1 

M 1-M. T 

(N.'1-.5)  F^An  (5,15'')  ((X  (I)  ,Y  (T)  ) ,L=5,M) 

PFAH  (5,15')  FB,  r.T,  SR  ,5ri,  XX,  YY  , XP  ,y  ,VY  , VT',  n'^I,  T>  ■ 

WRTTF  (5,250)  X?,YY',XX,YY,r.'5,5?,5F,E'',»^R,V~  r>i;'  ' 

WRITF  (5,45'')  ( (X  (T)  , Y (T)  ) ,T=1  , V)  ' 

no  5 t=1,25 
5 TYIJ  FT  (I)  ='*• 

nc  50  IOOP=1,2'''' 

XA=XX 
YA=YY 
TK  = 5 

C TNT'^TAl  T7.E  ALI'/ED  2AKDOF  Mn'^aFR,  0’-Nr;<A''0H 

r 

CALL  TIA  ’2  (7,'' ,rMT'^> 

TNT'"=IVTT*2+' 1 
CAT!.  CFAMD1  (TVIT) 

f 

r GFT  PRC  OIL 

DC  in  t-i,n 

(T)  ^n(X  (I)  ,Y  (I)  ,XF,YP)  /V2 
XF(T)  =X  (I)  fOPANn(SS,n.) 

YF  (3)  =Y  (T)  +GRANn(Sr.,*'  .) 
ir  CONTTNUF 

r 

C RFTMRfl  ’='TPF  ""O  (VR,YA) 

C 

5''  TK=TK  + 1 

IF  (TK.GT.25)  GOTO  "0 
TF=  " 

<T=XA  +001 
YT=v  A + DFT, 

^?  = FPAND  (5P, 

Y^D=GP  AMD  (Sn,r  .) 

■>MT-A  TAN  (XT/YT) 
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pr.T-Pi/0.  -rqp 
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- TP'"  riirnrT  '"pavpi.  tiffs 
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on  I'  vi-i,N 

'’T  (J)  =P(Xn)  ,Y  (.T)  ,Xr,vj)/V7 
SF(J)  =(YT-YF  (.1)  )/(XT-VF  (.1)) 

(SM  (J)  ) 

7F  ('T'H.L  T. ) Tll  = T'lfr>T 
r (J)  =C05  (TM) 

■'  (J)  =S’N  (T7) 

S'M.1)  =-i  ./0^‘  (‘^) 

V 70NTIKir 

(1)  (1)  -"-R  (1) 

K (?)  ='^T  { >)  (?) 

PO  P-1,’'' 

70  2''  L=1,7 

''=TT  (’■)  -'T  (T)  ♦P'lF i*qp A’cn  r .)  /’  ♦"  (f.i 

.7"  ' (P  , L ) =y  T-  "P*  VF*  s (I)  (T)  ”1-  S"  (■  ) ■"*VF  •i'C  ( P ) 

or  6 •'  T = 4 , 
on  7 - T - 1 , 

XI  (T.)  =(F<  (T,I.)  ,I  n / (.IF  (T  - 1 ) -S'M”)  ) 

y 1 (T.)  =SF  (T-1)  *x  1 (L)  +p  (1-1  ,T  ) 

< :Mt  ) = CMT  ^ r M -B  (P  , T ) ) / ( S'^  ( P)  - r.v  (I  f 1 ) ) 

Y2  (L)  =F'*  (T)  *X2  (L)  +«  (T  ,L) 

7"  CnN""TNnF 

PF  (X  1 (?)  .FO.  Y1  (1)  ) ''^OTO  RP 
oFI  = (Y1  (7)  - Y1  (1)  )/  (X  1 (7)  - XI  ( 1)  ) 

TF  (K2(7)  ."’0.X2  (1)  ) S0""o  0" 

"HPT=  (Y2  (2)  -Y  2 (1)  ) /(X  2 (?)  -x  : (1)  ) 

Y'CFF  = (Y1  (1)  - .SMT^'XI  (1)  - V 2 (1)  ♦ .SFTP*y2  (1)  ) /(:>F  p*  I) 
''Tpf*=s’'T*y'’’"i-*-yi  (1)  -Fiy^xi  (i) 
snin  ss 
qr  <TP''-X1(?) 

'P  (X2  (7)  . "■Q.  X ’ (1)  ) FnTO  qp’ 

SPTT  = (Y2  (2)  - Y?  (1)  ) / (X2  (2)  -Xi  (1)  ) 

Y-rv=.svTT*yf’P’i4.Y2  (1)  - p,'-TT*y  2 (1) 

’ o O r 

t'  VTrrM-X  ('') 

PP  (X  1 (.' ) . ■' 1 . X 1 p ) ) GRPn  qp' 

'1YT  = (V’  (7|  - vl(1))/(Xl  (7)-Xl  (1)) 
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rOWMlRLt  WILLOW  RUN  LAaORATOMItS,  THE  UNIVinaiTV  OF  MICHIOAN 


TABLE  B.4 

SENSOR  COORDINATES 
SENSOR  POSITIONS  1 THROUGH  9 
(x.y) 

CAMP  PENDLETON 

(3906»8862) 

(3911,8609) 

(3817,8380) 

(3946,7905) 

(3967,7656) 

(3909,7386) 

(3882,7037) 

(3831,6793) 
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ronMIRLV  WILLOW  RUN  LABOMATUmKM.  TMf  UKiVKB»(TV  OF  MtC.HlOAN 


Seismometer  position  B 


Seismometer  position  7 


Seismometer  position  5 


Seismometer  position  4 


Seismometer  position  3 


Seismometer  position  2 


FIGURE  B.l.  VERTICAL  AXIS  SEISMOGRAMS  FROM  EXPLOSION  EVENT  23 
TWENTY  NINE  PALMS,  1.5  km  Site,  Scale  10  in/sec 
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Seismometer  position  8 


Seismometer  position  7 


fONMIRLI  WILLOW  RUN  LABORATOBIIS,  1 HI  UNIVKRSITY  OF  MICHIGAN 


Seismometer  position  6 


J V 


Seismometer  position  3 

— /'X.  y' 


Seismometer  position  2 


FIGURE  B.2.  VERTICAL  AXIS  SEISMOGRAMS  FROM  EXPLOSION  EVENT  24, 
TWENTY  NINE  PALMS,  1.5  km  Site,  Scale  10  in/sec 
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FQMMCALV  WILLOW  RUN  LABORATORIEA.  THE  UNIVERSITY  OF  MICHIGAN 


Seismometer  position  8 


Seismometer  position  7 


Seismometer  position  6 


Seismometer  position  5 


Seismometer  position  A 


Seismometer  position  3 




Seismometer  position  2 


FIGURE  B.3.  VERTICAL  AXIS  SEISMOGRAMS  FROM  EXPLOSION  EVENT  69, 
TWENTY  NINE  PALMS,  1.5  km  Site,  Scale  10  In/sec 


“oRMIrTTwiLLOW  run  l.*BO«ATOHIt8,  THE  UNIVtRSITV  OF  MICHIGAN 


Seismometer  position  2 


Seismometer  position 


Seismometer  position  4 


Seismometer  position  5 


Seismometer  position  6 


Seismometer  position  8 


FIGURE  B.4. 


VERTICAL  AXIS  SEISMOGRAMS  FROM  EXPLOSION  EVENT  71 
TWENTY  NINE  PALMS,  1.5  km  Site,  Scale  10  in/sec 
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FOMtMCALV  WILLOW  RUN  L.A»ORATOMlC».  TMt  UNIVK.7MTY  OK  MICHIGAN 


Seismometer  position  7 


FIGURE 


VERTICAL  AXIS  SEISMOGRAMS  FROM  EXPLOSION  EVENT  16 
TWENTY  NINE  PALMS,  1.0  km  Site,  Scale  2 in/sec 
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Seismometer  position  7 


Seismometer  position  6 


Seismometer  position  5 


Seismometer  position  4 


Seismometer 
position  3 


' ronMlRLYWIULOW  BUNLABORATORim.  THt  UNIVCRSITY  OF  MICMIOAN 


FIGURE  VERTICAL  AXIS  SEISMOGRAMS  FROM  EXPLOSION  EVENT  21, 

TWENTY  NINE  PALMS,  1.0  km  Site,  Scale  2 in/sec 
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Seismometer  position  7 


Seismometer  position  6 


Seismometer  position  5 


Seismometer  position  4 


Seismometer  position  3 


FIGURE  B.IO.  CORRELATION  FUNCTIONS,  EVENT  11  VERSUS  EVENT  10, 
TWENTY  NINE  PALMS,  1.0  km  Site,  Scale  15  in/sec. 
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ronMIWLV  WILLOW  PUN  LABONATORIIl.  THC  UNIVtPlITr  OF  UICHIOAN 


Seismometer  position  4 


FIGURE  B.ll.  CORRELATION  FUNCTIONS,  EVENT  11  VERSUS  EVENT  16, 
TWENTY  NINE  PALMS,  1.0  km  Site,  Scale.  15  in/sec 


rOKMEHLV  WILLOVK  BUN  LABOBATOBIt*.  THE  UNIVSBBITY  OB  MICMIOAN 


Seismometer  position  7 


FIGURE  B.12. 


COPRELATION  FUNCTIONS,  EVENT  11  VERSUS  EVENT  21, 
TWENTY  NINE  FLAMS,  1.0  km  Site,  Scale  15  in/ sec 
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POUMtflLV  WIU.OW  LAHOAATOniCn,  THK  UNIVCRtlTV  f3f  Mir.HIQAM 


Seismometer  position  9 


Seismometer  posltion*6  .a  a ^ k . 




Seismometer  position  8 


PORMCIILV  WILLOW  RUN  LAKMATORKS.  THE  UNIVEKBITlf  Or  MICHIOAR 


Seismometer  position  9 
Seismometer  position  8 

k 


If  V I ’ I V V 

Seismometer  position  7 


Seismometer  position  6 


Seismometer  position  5 


Seismometer  position  A 


Seismometer  position  3 


Seismometer  position  1 


FIGUl 


VERTICAL  AXIS  SEISMOGRAMS  FR(M  MORTAR  RECOIL  EVENT  155, 
CAMP  PENDLETON,  SCALE  iO  in/ sec 


Seismometer  position  9 


Selsmoweter  posrtMn  8 


^elsmometlCT  position 


Seismometer  position  6i 


\yv 


Seismometer  position  5 


Seismometer  position  4 


Seismometer  position  3 


Seismometer  position  2 


Seismometer  position  ]. 
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FIGURE  B.17.  AUTO-CORRELATION  FUNCTIONS,  EXPLOSION  EVENT  103, 
CAMP  PENDLETON,  SCALE  15  in/sec 


Seismometer  position  9 

\[\f\jA/\/'^^\^^^ 


Seismometer  position  6 


CAMP  PENDLETON,  SCALE  15  in/ sec 


Seismometer  position  9 


Seismometer  position  8 


Seismometer  position  6 


P‘, 

\i 

w 


1 


Seismometer  position  3 


CAMP  PENDnETON,  SCALE  15  in/sec 
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